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Program  Results 

Executive  Summary: 

4)  Statement  of  the  problem  studied 

The  goals  of  this  program  were  to: 

1 .  Provide  proof-of-principle  for  the  feasibility  of  using  soliton  collisions  for 
computing 

2.  Investigate  new  fundamental  concepts  in  soliton  science  to  address  #1. 

3.  Based  on  new  science  from  #2,  investigate  new  possible  applications  of  solitons. 

4.  Investigate  new  related  ideas. 

5)  Summary  of  the  most  import  results 

1 .  Soliton  Computing 

A  number  of  logic  functions  and  mathematical  operations  were  implemented  in 
the  laboratory  based  on  soliton  collisions  in  photoreffactive  media.  In  addition  to  the 
usual  NAND  and  AND  logic  gates,  soliton  collisions  do  transfer  information  and  two 
successive  collisions  can  be  made  to  mimic  a  unitary  matrix  or  its  inverse. 

Proof-of-principle  investigations  demonstrated  that  Manakov  vector  soliton-based 
computing  is  Turing  Universal:  any  computation  can  be  implemented  in  a  system  of 
colliding  solitons.  Furthermore,  because  these  devices  are  based  on  a  reversible  physics, 
such  a  computing  system  would  not  generate  waste  heat  in  the  collisions  that  implement 
its  bit-level  operations.  Our  comparative  analysis  pre-supposed  a  performance  target  of 
one  Peta  FLOPS,  and  compared  the  efficiency  of  implementations  in  solitonic  and 
standard  silicon  VLSI  technologies.  While  we  have  demonstrated  a  clear  proof-of- 
principle  for  soliton  based  optical  computing  functions,  electronics  has  advanced  in  speed 
and  power  requirements  in  the  last  six  years  and  we  recognize  that  there  are  power 
requirements  that  still  need  to  be  addressed  in  the  solitonic  case. 

2.  A  plethora  of  new  solitons  and  soliton  phenomena  were  discovered  in  the  course 
of  this  program.  In  continuous  and  discrete  media,  the  concept  that  solitons  need 
coherence  was  shattered  and  a  variety  of  new  solitons  were  reported  including  random 
phase  solitons.  Solitons  with  unique  properties  were  demonstrated  in  semiconductor 
optical  amplifiers  at  lOmW  power  levels.  But  by  far  the  most  interesting  results  were 
obtained  in  discrete  systems  both  for  single  peaked  and  vortex  solitons.  Because  the 
diffraction  in  such  systems  can  be  either  normal  or  anomalous  many  new  types  of 
solitons  were  obtained  in  photorefractive,  Kerr  and  quadratic  nonlinear  media.  In  addition, 
solitons  in  multiple  bands  were  observed  and  collided.  Parallel  studies  of  filamentation 
indicated  that  there  were  regions  of  operation  in  which  filamentation  did  not  occur. 
Another  highlight  was  the  Brillouin  zone  spectroscopy  by  which  the  Fourier 
representation  of  the  Brillouin  zone  can  be  mapped  out  on  a  single  experiment.  The 
richest  results  by  far  occurred  in  this  basic  science  part  of  the  program. 

3.  A  number  of  new  device  concepts  emerged  from  the  basic  soliton  physics  studies 
in  2. 
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(1)  Since  the  output  of  the  discrete  arrays  is  spatially  digital,  a  digital  beam  scanner 
was  demonstrated. 

(2)  Furthermore,  it  was  proposed  that  an  all-optical  2D  router  could  be  implemented 
based  on  collisions  between  highly  localized  “blocker”  solitons  and  beams  traveling  in 
diffraction  free  directions,  or  weakly  localized  solitons.  The  basic  concepts  were  verified 
in  ID  arrays  with  the  incoherent  interaction  being  the  most  attractive. 

(3)  One  of  the  most  interesting  outcomes  has  been  the  development  of  an 
understanding  of  nano  optical  waveguides  and  their  potential  benefits  to  Army  needs. 

4.  Fabrication  of  an  important  nonlinear  material,  namely  quantum  dots,  was 
initiated  and  significant  progress  made.  Key  results  such  as  control  of  quantum  dot 
topology  were  achieved.  In  a  separate  development,  digital  operation  of  SOAs  was 
demonstrated  as  triggerable  memory  elements. 

This  program  dictated  and  led  developments  in  the  field  of  soliton  science  and  its 
applications.  Its  scientific  impact  was  immense.  It  can  be  quantified  in  terms  of  the 
publications  and  presentations  that  flowed  from  it.  Specifically,  98  of  the  total  of  230 
published  papers  in  refereed  journals  were  in  the  high  impact  factor  letters  journals 
Science,  Nature,  Physical  Review  Letters,  Optics  Letters,  Optics  Express  and  Applied 
Physics  Letters.  Finally  there  were  98  invited  papers  at  international  conferences,  of 
which  4  were  plenary.  In  addition  there  were  187  contributed  conference  presentations. 

Under  sponsorship  of  this  MURI  program  there  were  16  PhDs  awarded  who  are 
now  in  the  workforce,  4  MScs  and  1 1  undergraduates  were  introduced  to  research. 
Finally,  there  were  a  total  of  1 1  senior  personnel  which  include  post-doctoral  fellows, 
research  associates,  visiting  faculty  etc. 
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Squier  -  Summary  Report  on  Soliton  Computing 

4. 1  Logic  and  General  Computation 

This  project  began  by  asking  the  basic  question,  Can  colliding  solitons  be  used  to  implement 
general  computation?  The  work  done  on  colliding  Manakov  vector  solitons  developed  a 
collection  of  digital  devices  (NAND  gates,  fanout,  and  signal  crossover)  that  taken  together 
form  a  Boolean  complete  set.  That  is,  using  only  these  devices,  one  can  implement  any 
arbitrary  Boolean  circuit.  These  devices  include  logic  level  restoration  and  clocked  state 
elements,  thereby  allowing  unbounded  cascading  of  devices  and  synchronous  system  design. 
With  these  elements  in  hand,  we  can  state  that  Manakov  vector  soliton-based  computing  is 
Turing  Universal:  any  computation  can  be  implemented  in  a  system  of  colliding  solitons. 
Indeed,  any  conventional  computing  machinery  organization  can  be  implemented  this  way. 
Furthermore,  because  these  devices  are  based  on  a  reversible  physics,  such  a  computing 
system  would  not  generate  waste  heat  in  the  collisions  that  implement  its  bit-level  operations. 

These  results  established  that  systems  of  colliding  optical  solitons  have  no  inherent  limitation 
restricting  their  ability  to  implement  computing  machinery.  However,  can  colliding  solitons 
be  used  efficiently  to  implement  computation?  For  a  first  approximation,  efficiency  is  in 
terms  of  processing  speed,  power  consumption,  and  material  resources  measured  in  area  of 
nonlinear  optical  propagation  media.  Of  course,  optimal  efficiency  depends  on  the  degree  to 
which  a  computation  makes  optimal  use  of  the  state  transformations  embodied  in  soliton 
collisions.  It  may  be  that  the  digital  logic  devices  described  above  are  substantially  sub- 
optimal,  and  some  other  way  of  employing  soliton  collisions  is  more  effective.  Naturally,  if 
some  application  algorithm  were  known  to  optimally  match  the  characteristics  of  a  system  of 
colliding  solitons,  implementing  that  algorithm  would  yield  the  system’s  efficiency. 

But  even  without  such  a  benchmark,  one  can  bound  computational  effectiveness  by  using  a 
comparative  analysis  that  looks  in  parallel  at  another  technology.  Our  comparative  method 
pre-supposes  a  performance  target  of  one  Peta  FLOPS,  and  compares  the  efficiency  of 
implementations  in  solitonic  and  standard  silicon  VLSI  technologies.  Similar  simplifying 
assumptions  are  used  for  both  implementations,  and  the  power  and  size  requirements  are 
compared.  Since  we  do  not  know  the  most  effective  method  of  using  soliton  collisions,  we 
look  at  a  collision  system  without  specifying  the  devices  or  system  structure. 

Estimates  of  the  performance  of  the  solitonic  system’s  internal  devices  are  based  on  the 
assumed  computational  complexity  of  a  single  collision.  In  the  least  restrictive  model,  we 
simply  count  the  maximum  number  of  collisions  possible,  and  use  an  absolute  upper  bound 
on  the  complexity  of  the  computational  operation  taking  place  in  a  single  collision.  By 
adding  more  restrictive  and  realistic  assumptions,  we  progress  from  model  to  model 
establishing  progressively  more  conservative  performance  bounds.  The  most  restrictive 
model  uses  a  simplification  of  an  implementation  of  a  floating-point  unit  in  digital  logic, 
which  is  replicated  as  needed  to  reach  the  performance  target.  Throughout  the  analysis,  the 
precise  nature  of  the  computation  is  left  unspecified,  but  it  is  assumed  there  is  always 
sufficient  instruction-level  parallelism  available  to  sustain  any  computational  units  at  full 
speed.  There  is  no  accommodation  for  such  things  as  pipeline  registers,  branching  hazards, 
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data  dependencies,  or  register  conflicts,  and  there  are  no  resources  devoted  to  system  control 
functions  or  instruction  flow. 

In  general,  we  consider  a  collision  system  with  n  spatial  2-component  Manakov  vector 
solitons  entering  and  exiting  a  collision  system  (see  Figure  1).  The  collisions  internally 
implement  a  computation  with  the  data  inputs  encoded  in  the  states  of  the  input  solitons  and 
the  system’s  output  encoded  in  the  states  of  the  output  solitons.  Not  every  input  soliton 
necessarily  carries  input  data,  some  may  be  needed  to  supply  non-data  inputs  to 
computational  collision  devices.  Likewise,  not  every  output  soliton  carries  output 
information,  some  may  carry  away  garbage  information.  But,  all  computation  is  internal  to 
the  optical  medium  and  is  in  the  form  of  soliton  collisions,  and  there  are  no  soliton  sources 
internal  to  the  medium. 


Figure  1 

A  collision  system  with  n  input  solitons  and  n  output  solitons.  There  are  no  internal  soliton  sources.  All 
computation  is  implemented  within  the  soliton  collisions  inside  the  optical  medium. 

The  maximum  number  of  collisions  occur  if  all  the  solitons  are  co-planar.  Each  soliton  is 
one  micron  wide,  with  two  micron  separation  between  non-colliding  solitons.  For  a 
typical  collision  angle  of  about  1°,  the  beam  crossing  distance  is  about  100  um  and  thus 
the  time  for  one  soliton  to  cross  another  soliton’s  path  is  about  0.5  ps  (see  Figure  2). 
Assuming  a  minimum  separation  between  collisions  of  two  beam  widths,  the  distance 
between  collisions  is  then  three  times  the  collision  length,  or  about  300  um.  The  total  area 
assigned  to  a  collision  accounts  for  a  2-um  separation  distance  between  solitons  paths.  A 
collision’s  delay  is  then  just  the  propagation  time  across  this  minimum  separation  plus  the 
beam  crossing  time.  Layout  and  signal  routing  overhead  is  accounted  for  by  using 
efficiency  factors  for  VLSI.  This  and  other  modeling  parameters  come  from  the 
International  Technology  Roadmap  for  Semiconductors  (1TRS).  Characteristics  for  our 
VLSI  model  are  estimated  using  the  1TRS  hp90  node,  which  is  current  technology  with  a 
feature  size  of  90nm. 
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Figure  2 

Three  colliding  solitons.  The  red  area  is  a  collision,  the  green  areas  are  associated  with  that  collision.  Total 
area  charged  to  each  collision  is  three  times  the  sum  of  the  red  and  green  areas.  Delay  time  charged  to  the 
collision  is  the  propagation  time  from  collision  to  collision,  which  is  three  times  the  soliton  crossing  time. 

Using  these  basic  parameters,  we  implemented  a  model  of  a  64-bit  adder/multiplier  unit  in 
VLSI  and  then  converted  to  equivalents  in  the  digital  solitonic  model  (see  Figure  3).  For 
instance,  while  the  VLSI  design  was  not  logic  gate  based,  we  used  the  ITRS  mean  of  four 
transistors  per  logic  gate,  and  converted  the  VLSI  transistor  count  to  an  equivalent  number 
of  solitonic  logic  gates.  The  design  was  fully  pipelined  in  a  manner  that  is  not  practically 
possible,  but  conforms  to  modeling  assumptions  made  for  optical  logic  in  the  literature. 


four  64-bit  inputs 


C  four  64-bit  outpus  ) 

Figure  3 

A  pipelined,  parallel,  floating-point  processing  unit.  Units  are  replicated  as  needed  to  attain  the  performance 
goal. 

That  is,  without  using  pipeline  registers,  the  unit  can  receive  new  data  inputs  after  waiting 
only  a  nominal  delay.  The  delay  for  VLSI  model  is  based  on  the  delay  of  a  fanout-of-four 
(F04)  inverter.  Six  F04  delays  is  chosen  as  that  is  an  apparent  limit  on  delay  for  the 
fastest  pipeline  design  possible.  Being  more  generous  in  the  solitonic  model,  we  use  the 
delay  of  two  basic  pseudo-logic  gates,  based  on  the  fact  that  any  logic  function  can  be 
implemented  in  two-layer  logic  if  one  ignores  fanout.  These  gates,  which  we  call 
Steiglitz-Rand  gates  (SR  gates),  are  a  composite  of  Steiglitz’s  NAND  gate  and  Rand’ 


5 


NAND  gate,  using  the  least  restrictive  assumptions  in  both  cases.  That  is,  we  combine  the 
best  performance  and  efficiency  features  from  known,  solitonic  digital  logic  devices 
without  regard  for  the  existence  of  any  logic  device  having  such  a  mix  of  characteristics. 
While  not  an  actual  logic  gate,  this  pseudo-logic  gate  models  a  cascadable,  signal¬ 
restoring  logic  with  the  least  delay  and  fewest  non-data  inputs.  In  general,  only  the  raw 
computing  power  based  on  device  count  is  considered,  and  only  optimistic  assumptions 
for  solitonic  computing  systems  are  used. 


Technology 

VLSI 

Solitonic-digital  1 

Solitonic-digital  2 

Feature  Size  (um) 

90  nm 

1000  nm 

100  nm 

Soliton  power 

10  mW 

25  uW 

Total  Area 

30  cm2 

1  m2 

23  cm2 

Total  Power 

1  kW 

5  MW 

1.5  kW 

Table  1 

Comparison  of  silicon-based  versus  solitonic-based  1-Peta  FLOPS,  64-bit  floating  point  compute  engines. 

The  results  are  shown  in  Table  1  for  implementations  of  a  1-Peta  FLOPS  compute  engine 
based  on  these  adder/multipliers.  The  sources  for  solitons  that  are  not  cascaded  between 
logic  gates  in  the  “Solitonic-digital  1”  model  consume  10  mW  each  and  produce  1-um 
solitons.  To  reach  the  breakeven  point  with  VLSI,  the  “Solitonic-digital  2”  model  has  25 
uW  soliton  sources  producing  100  nm  solitons. 

It  is  important  to  note  what  this  analysis  does  and  does  not  do.  We  do  not  start  with 
optical  Boolean  logic  devices,  build  a  computing  machine,  and  then  measure  its 
performance.  Rather,  we  first  design  a  hypothetical  model  of  a  VLSI  system  that 
consumes  the  least  resources  and  has  the  best  performance.  The  model  serves  to  identify 
the  minimum  resources  required  for  a  system  with  the  target  performance.  Simplifying 
assumptions  are  always  in  favor  of  optimistic  efficiency  and  speed  over  realistic 
possibilities.  This  gives  a  lower  bound  on  required  resources  with  an  upper  bound  on 
performance.  The  equivalent  optical  model  is  then  constructed  using  an  optimistic  bit- 
operations-per-soliton-collision  factor.  This  conversion  factor  is  derived  for  the  digital 
solitonic  model  using  the  SR  gate.  Both  models  then  are  used  to  produce  area  and  power 
requirements  based  on  current  practice  for  VLSI  and  hypothetical  bounds  for  solitonic 
systems. 

The  result  is  a  picture  of  what  kind  of  resources  are  required  to  achieve  a  given  level  of 
performance.  The  picture  is  grossly  over-optimistic.  But,  it  approximates  real  world 
expectations  to  the  degree  that  actual  system  design,  efficiency,  and  performance 
parameters  have  been  preserved  in  the  modeling.  So,  it  presents  a  picture  of  solitonic 
performance  and  efficiency  that  is  relatively  easy  to  understand  since  it  gives  parameters 
of  a  vaguely  realistic  computing  machine.  The  comparative  analysis  also  provides  a 
check  on  the  mismatch  between  reality  and  the  models.  Clearly,  no  one  would  expect  to 
produce  a  1-Peta  FLOPS  adder/multiplier  using  only  ten  300  mm2  silicon  chips. 

However,  the  comparison  between  the  two  technologies  shows  what  kind  of  solitonic 
power  and  area  efficiencies  would  be  needed  to  make  solitonic  systems  competitive  for 
computing,  under  the  assumption  of  binary  digital  processing. 
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Salamo  -  Arkansas  -  Summary  Report 
I.  Solitons  for  Computing 

A.  Basic  Science:  During  this  MURI  we  experimentally 
demonstrated  second-band  modulation  instability  (MI)  in  a  self- 
focusing  waveguide  array,  obtained  by  exciting  modes  at  the  edge  of 
the  first  Brillouin  zone.  This  observation  was  directly  compared  with 
first  band  MI  at  the  zone  edge  by  switching  to  a  defocusing 
nonlinearity.  In  both  cases,  the  nonlinearity  pushes  the  propagation 
constant  off  the  linear  bands,  i.e.  these  nonlinear  modes  reside  within  a 
band  gap  of  the  transmission  spectrum.  That  is,  we  have  demonstrated 
lattice  MI  in  photonic  band  gaps.  Moreover,  the  observations  are 
unique  to  periodic  systems  and  do  not  occur  in  homogeneous  media. 

We  also  observed  second-band  MI,  the  corresponding  underlying 
theory,  and  distinguish  the  observed  features  from  first-band  MI  and 
MI  in  homogeneous  media.  For  the  results  in  Figure  1,  with  intensity 
ratios  I  =  0.1  and  0.25,  the  experimental  MI  period  corresponds  to 
spatial  momentum  0  =  0.56  and  0  =  0.78,  respectively,  in  good 
agreement  with  the  corresponding  theoretical  values  0  =  0.52  and  0  = 

0.80.  We  also  note  that,  similar  to  band  1,  the  MI  gain  and  spatial 
frequency  for  band  2  increase  with  higher  intensity. 

B.  Applied  Science:  One  of 
the  principal  goals  of  this  MURI  was 
to  use  soliton  collisions  in  an  efficient 
way  to  do  computation  in  a  flexible 
manner.  It  was  clear,  theoretically  and 
experimentally,  that  soliton  collisions 
were  a  viable  way  to  transfer 
information  using  vector  solitons. 

A  vector  soliton  is  a  light 
beam  composed  of  two  independent 
components  (A,B).  During  this  MURI 
multiple  collisions  between  vector 
solitons  can  be  used  to  develop  a 
Matrix  Algebra. 

In  particular,  we  have 
demonstrated  that  two  vector  solitons 
can  act  as  either  a  unity  matrix  that 
gives  back  the  same  vector  after  the 
matrix  two  collisions  or  the  inverse 
matrix  that  reverse  the  A  and  B 
components.  In  the  figure  below  we 
show  that  the  incoming  vector  is 
made  of  A  and  B  components  of  ratio 
Al/Bl  =  0.49  and  that  the  Matrix 
operation  reverses  the  two 
components  so  that  the  ratio  Bl/Al  = 

0.49.  Of  course,  these  operations 
conserve  the  magnitude  of  each 
vector  so  that  A1  +  B1  remains 
constant. 

Fig.  1.  (a)  Unity;  (b)  inverse  operations 
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FIG..  Observation  of  lattice 
MI  for  a  defocusing  ID  array 
for  (a)  1=0.1  and  (b)  I  =  0.25 , 
and  (c)  when  the  F-B  mode 
located  at  the  edge  of  band-2  is 
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II.  Quantum-Dot  3D  Lattices  for  Enhanced  Nonlinear  Optical  Soliton  Interactions 

A.  Basic  Science:  The  quest  for  computing  concepts  requires  new  approaches  to  material  growth 
as  well  as  new  materials  systems.  During  this  MURI  we  have  developed  proficiency  in  molecular  beam 
epitaxy  (MBE),  electron  microscopy,  and  optical  probe  techniques  which  resulted  in  our  published 
success  in  the  fabrication  and  understanding  of  the  behavior  of  self-assembled  nano  structures.  Building 
on  this  early  foundation  during  the  first  two  years  of  the  MURI,  we  have  been  advancing  the  science  and 
tailoring  the  behavior  of  organized  arrays  of  nanostructures.  These  advances  provide  the  basis  for 
crafting  new  photonic  lattices  with  enhanced  nonlinear  behavior.  Our  goal  of  the  MURI  was  to  refine  our 
skills  in  MBE  to  yield  beautifully  ordered  2D  and  3D  arrays  of  quantum  dots,  wires,  and  rings  (noted  in 
Science).  Achieving  this  control  over  growth  is  now  yielding  material  that  gives  new  insight  into  the 
collective  interactions  between  individual  units  and  provides  the  basis  to  tailor  a  number  of  remarkable 
properties  of  3D  arrays:  geometry  dependent  excited  state  lifetimes;  tailored  refraction  and  dispersion, 
and  enhanced  nonlinear  optical  coefficients. 

As  mentioned  during  the  first  two  years  of  the  MURI  we  used  our  MBE  and  in-situ  STM  and 
explored  the  nature  of  the  size,  shape,  density,  aspect  ratio  and  volume  of  III- V  semiconductor  3-D  self- 
assembled  islands.  For  example,  our  data  for  3-D  InAs  islands  on  a  GaAs  substrate  revealed  for  the  first 
time  the  structural  details  of  two  different  shaped  islands.  These  pictures  (Fig.  1)  beautifully  demonstrate 
the  power  of  STM  to  reveal  shape  and  size  of  quantum  dot  structures. 


Fig.  1.  STM  picture  and  structure  for  two  different  dot  shapes  observed 


More  recently,  however,  we  have  explored  their  collective  behavior.  Periodic  nanostructure 
arrays  are  desired  for  two  reasons.  First,  they  are  much  more  uniform,  which  greatly  improves  their 
efficiency  as  lasers,  detectors,  memory  systems,  and  nonlinear  optical 
devices.  Second,  the  collective  properties  of  these  arrays,  such  as  their 
index  of  refraction,  their  polarizability,  or  their  band  structure  depend  upon 
their  spacing,  degree  of  interaction,  defects,  characteristic  single  electron 
energies,  etc.  While  we  are  well  aware  of  the  collective  behavior  of  solids, 
we  have  a  new  playing  field  to  study  and  engineer  “collective  behavior”. 

Shown  in  Fig.  2  is  the  x-ray  diffraction  pattern  of  2D  lattices  of  QDs.  The 
peaks  in  the  x-ray  pattern  tell  us  that  we  have  an  array  that  mirrors  a  typical 
crystal.  As  a  result,  we  can  expect  that  by  nano-managing  the  structure  of 
an  organized  array  of  synthetic  units  we  can  create  systems  in  which  we 
have  control  of  their  optical  and  electronic  properties,  even  tailoring  them 
to  regimes  not  found  in  nature.  For  example,  we  are  creating  media  with 
abnormally  large  or  even  negative  indices  of  refraction  or  nonlinear  optical 
photonic  lattices  that  can  dramatically  change  their  reflectivity  with  an 
applied  electric  field  or  increasing  optical  intensity.  Such  “material  by 
design”  control  would  open  a  new  era  in  applications  ranging  from 


Fig.  2.  An  organized  QD 
2D  array.  The  insert  is  an  x- 
ray  pattern  that  reveals  the 
array  periodicity.  QD  size  is 
20  to  30  nnt. 


improving  nonlinear  coefficients  to  optical  computation  that  rival  their  electronic  counterparts. 
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Fig  3  (a).  ID  array  of 
GaAs/InGaAs  QDs 


Fig.  3(b).  Cross-section  of  a  quantum  dot  3D  array 


In  addition  to  2D  quantum  dot  nanostructures  we  have  developed 
beautiful  long  chains  of  quantum  dots  that  were  highlighted  in  Science. 

These  chains  are  amazingly  over  five  microns  long  and  have  improved 
homogeneity  when  compared  to  the  growth  of  randomly  order  dots  (Fig.  3 
(a)).  These  arrays  already  show  interesting  organizational  as  well  as 
optical  and  transport  behavior  and  offer  the  interesting  opportunity  to 
explore  dots  chains  as  electrical  nano-interconnects.  They  are  fabricated 
by  layering  quantum  dots  as  seen  in  the  TEM  picture  below  (Fig.  3  (b)). 

The  picture  gives  atomic  resolution  so  that  we  can  also  measure  the 
quantum  dot  facets  (see  Phy  Rev  B  69,233312  (2004);  APL  85,  5031 
(2004);  APL  84,  1931  (2004);  APL  84,  4681  (2004)).Using  our 
experience  in  coherent  pulse  propagation  and  quantum  optics  to  design 
and  measure  dot  arrays  with  tailored  absorption  and  dispersion  properties 
and  even  nonlinear  coefficients. 

Moreover,  using  electromagnetic 
induced  transparency  (EIT) 
technology  we  will  be  able  to 
selectively  excite  different  regions 
or  clusters  of  nanostructures 
creating  dynamic  interference  effects  that  can  be  used  to  control  the  magnitude  of  the  optical  nonlinearity. 

In  this  MURI  we  have  also  worked  on  semiconductor  ring 
structures  as  small  as  40nm  across  will  be  fabricated  via  self-assembly 
(Fig.  3  (c)).  Ring  structures  differ  from  dots  in  that  a  classical  point 
charge  on  a  ring  always  has  a  dipole  moment  of  fixed  magnitude, 
although  of  arbitrary  direction.  Theory  indicates  that  nano  ring  arrays 
should  therefore  undergo  a  phase  transition  from  isotropic  to  anti- 
ferroelectric  ordering  as  their  interactions  are  increased.  The  point  here  is 
that  it  may  be  possible  to  tune  ring  geometries  and  produce  structures 
with  anisotropic  polarizabilities  and  larger  nonlinear  coefficients.  Such 
materials  may  have  applications  for  switching  optical  solitons. 

B.  Applied  Science:  Finally,  we  have  just  started  to  grow 
ferroelectric  nanostructures  ranging  from  slabs,  to  rods  to  dots.  The 
anticipated  benefit  hinges  on  whether  the  phase  transitions  and  multi¬ 
stable  states  still  exist  in  low-dimensional  structures.  This  question  is  the 
subject  of  much  contention  and  long-standing 
interest  for  understanding  and  revealing 
collective  interactions.  With  many  basic 
questions  to  answer  and  significant  potential 
application,  there  are  tantalizing  opportunities 
for  theory-experiment  to  push  the  state-of-the- 
art.  The  precise  effects  of  the  substrate,  growth 
orientation,  surface  termination  and  thickness 
on  the  ferroelectric  and  ferromagnetic 
properties  of  nanodot  arrays  are  also  not  clear. 

What  will  be  their  linear  and  nonlinear  optical 
character  and  can  they  be  enhanced?  We  have 
the  unique  opportunity  to  uncover  the  nature  of 
nanoferroelectrics  and  increase  their  potential 
for  application.  In  this  MURI  we  have  already 
reported  (Fig.4)  enhanced  nonlinearity  in  our 
nanostructure  Bao.6Sr0.4Ti03  grown  films. 


Fig.  3  (c)  Arrays  of 
GaAs/AlGaAs  QRs 
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Fig.  4  Z-scan  Measurements  of  Ba06Sr04TiO3  film  on 
MgO  at  the  wavelength  of 395  nm 
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Summary  Report:  Synopsis  of  MURI  activities 
Christodoulides-University  of  Central  Florida 

During  this  MURI,  my  group  has  been  engaged  in  theoretical  and  modeling  activities 
concerning  soliton  interactions.  In  particular  we  have  investigated  the  possibility  of 
utilizing  discrete  solitons  in  waveguide  arrays  and  lattices  for  all-optical  switching  and 
routing  applications.  This  activity  was  carried  out  in  very  close  collaboration  with  the 
experimental  groups  involved  in  this  MURI.  In  what  follows  we  highlight  some  of  the 
main  thrusts  of  these  activities. 

Solitons  for  signal  processing  applications: 

1.  Navigating  discrete  solitons  in  two-dimensional  nonlinear  waveguide  array 
networks: 

Among  the  most  important  functions  of  a  photonic  network  is  to  route  information  from  a 
particular  point  of  origin  A  to  its  final  destination  port  Z.  In  such  optical  systems  it  is 
often  highly  desirable  that  routing  is  accomplished  all-optically  so  as  to  avoid 
unnecessary  electro-optic  conversion.  If  for  example  data  is  re-directed  by  a  space¬ 
switching  matrix,  it  is  also  crucial  that  this  process  occurs  with  minimum  diffraction 
induced  cross-talk  or  losses  among  nodes.  Nonlinearity  can  offer  a  promising  solution  to 
this  problem  since,  under  appropriate  conditions,  is  known  to  balance  diffraction  effects. 
In  fact,  in  nonlinear  waveguide  arrays  (which  can  serve  as  the  network  nodes)  a  self- 
trapped  entity  is  possible,  better  known  as  a  discrete  soliton  (DS).  By  their  very  nature, 
discrete  solitons  represent  collective  excitations  of  the  periodic  lattice  as  a  whole  and 
they  are  the  outcome  of  the  balance  between  nonlinearity  and  linear  coupling  or  discrete 
diffraction  effects.  Lately,  optical  DS  were  successfully  observed  in  nonlinear  AlGaAs 
waveguide  arrays  (by  George  Stegeman  and  Yaron  Silberberg).  As  noted  in  several 
studies,  the  discrete  character  of  these  self-trapped  states  can  lead  to  a  host  of  new  effects 
that  have  no  counterpart  whatsoever  in  the  continuous/bulk  regime. 

In  a  series  of  papers,  we  have  shown  that  DS  in  two-dimensional  nonlinear 
waveguide  array  networks  can  provide  a  rich  environment  for  all-optical  data  processing 
applications  (PRL  87,  pp.  233901,  2001).  More  specifically  we  have  demonstrated  that 
this  family  of  solitons  can  be  employed  to  realize  intelligent  functional  operations  such 
as  routing,  blocking,  logic  functions,  time-gating  etc.  These  DS  can  be  navigated 
anywhere  in  the  network  along  pre-assigned  array  pathways  which  act  like  “soliton 
wires”.  Even  more  importantly,  DS  can  be  routed  at  array  intersections  using 
vector/incoherent  interactions  with  other  discrete  solitons.  In  essence,  these  intersections 
behave  as  DS  switching  junctions. 

Figure  1(a)  depicts  a  nonlinear  array  network  involving  consecutive  bends.  Light 
in  this  array  propagates  along  the  z-axis  and  is  confined  in  the  transverse  x-y  plane.  The 
waveguide  cross-sections  are  shown  in  green.  A  discrete  soliton  S,  the  field  profile  of 
which  is  shown  in  red,  is  set  in  motion  in  this  system.  This  is  done  by  appropriately 
chirping  (spatially)  or  tilting  the  soliton  beam  with  respect  to  the  z-axis.  In  order  to 
investigate  the  effects  of  bends  on  the  behavior  of  DS,  we  have  numerically  simulated  the 
process.  Interestingly  enough  we  have  found  that  DS  can  successfully  negotiate  a 
sequence  of  bends  with  very  little  radiation/reflection  losses.  The  DS  follows  the  pre- 


10 


assigned  path  (as  shown  schematically  in  Fig.  1(a))  and  remains  essentially  invariant 
during  propagation. 


Such  systems  can  also  be  used  to  realize  logic  operations.  Figure  1(b)  illustrates  an  X  - 
switching  junction.  This  junction  employs  two  different  DS  families,  the  so-called 
“signals”  and  “blockers”  denoted  by  S  and  B,  respectively.  Signals  are  moderately 
confined  DS  whereas  blockers  (depicted  in  blue)  are  strongly  confined  occupying 
effectively  one  site.  Unlike  signals  which  are  highly  mobile,  blockers,  states  deeply 
located  in  the  spatial  photonic  bandgap  of  the  waveguide  lattice,  tend  to  retain  their 
position  after  a  collision  event.  In  the  junction  shown,  the  blockers  Bl  and  B2  interact 
“incoherently”  with  S.  This  is  possible  by  using  different  colors  or  polarizations  in  a 
nonlinear  material.  As  the  figure  schematically  shows,  the  signal  soliton  S  is  routed 
towards  the  lower  branch,  due  to  the  presence  of  the  two  blockers  at  the  entries  of  the 
respective  pathways.  This  again  happens  with  very  little  loss  (below  4%  per  junction).  It 
is  important  to  note  that  had  one  of  the  two  blockers  not  been  present  at  the  junction,  the 
signal  DS  would  have  totally  disintegrated  into  transmitted  and  reflected  waves.  Thus  in 
essence,  the  junction  operates  as  an  AND  gate.  The  design  of  such  DS  switching 
junctions  was  reported  in  Optics  Letters,  vol.  26,  pp.  1978  (2001).  Similarly,  other 
operations,  such  as  time-gating  and  memory  functions,  can  also  be  implemented  in  these 
systems. 

2.  Minimizing  bending  losses  in  discrete  soliton  networks: 

Bends  in  discrete  soliton  networks  are  certainly  essential  elements  in  such  array  systems. 
In  general,  DSs  propagating  along  a  bend  tend  to  lose  power  because  of  reflections. 
Depending  on  the  bend  angle,  these  losses  can  be  considerable.  This  is  especially 
problematic  if  bends  are  to  be  used  in  cascade.  Reducing  bending  losses  in  DS  networks 
is  thus  an  issue  of  importance.  In  this  period  we  have  shown  (Optics  Letters,  vol.  26,  pp. 
1876)  that,  by  appropriate  design,  reflection  losses  occurring  along  sharp  bends  in  two 
dimensional  DS  networks  can  be  almost  eliminated.  Analysis  indicates  that  this  can  be 
accomplished  by  modification  of  the  guiding  properties  of  the  comer  waveguide.  Our 
analytical  results  were  found  to  be  in  excellent  agreement  with  numerical  simulations. 
The  design  of  these  comers  was  carried  out  by  employing  the  corresponding  unfolded 
photonic  circuit  of  the  bend.  Similarly,  the  switching  efficiency  at  T  and  X-junctions  can 
be  optimized  by  engineering  the  waveguides  at  the  vicinity  of  the  switching  junction  as 
shown  in  the  figure  below. 
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The  performance  of  switching  junctions  in  two-dimensional  discrete-soliton  networks 
was  analyzed  theoretically  by  coupled-mode  theory.  Our  analysis  can  be  used  for  the 
design  of  routing  junctions  with  specified  operational  characteristics.  Appropriately 
engineering  the  intersection  site  can  further  improve  the  switching  efficiency  of  these 
junctions  up  to  99%.  Our  analytical  results  were  verified  again  by  numerical  simulations. 

3.  Discrete  temporal  solitons  along  a  chain  of  nonlinear  coupled  microcavities 
embedded  in  photonic  crystals: 

We  have  shown  that  spatiotemporal  discrete  solitons  are  possible  in  nonlinear  photonic 
crystal  structures.  Analysis  indicates  that  these  states  can  propagate  undistorted  along  a 
series  of  coupled  resonators  or  defects  by  balancing  the  effects  of  discrete  lattice 
dispersion  with  material  nonlinearity.  In  principle,  these  self-localized  entities  are  capable 
of  exhibiting  very  low  velocities  depending  on  the  coupling  coefficient  among  successive 
microcavities.  This  class  of  solitons  can  follow  any  pre-assigned  path  in  a  three- 
dimensional  environment.  These  results  have  been  reported  in  Optics  Letters,  vol.  27,  pp. 
568  (2002). 

4.  Discrete  solitons  in  nonlinear  zigzag  optical  waveguide  arrays: 

We  have  demonstrated  that  the  discrete  diffraction  properties  of  a  nonlinear  zigzag 
waveguide  array  can  be  significantly  modified  by  exploiting  the  topological  arrangement 
of  the  lattice  itself.  This  introduces  long-range  interactions  which  in  turn  affect  the  lattice 
dispersion  relation  within  the  Brillouin  zone.  As  a  result  of  this  band  alteration  we  have 
demonstrated  that  new  families  of  discrete  solitons  are  possible  which  are  stable  over  a 
wide  range  of  parameters.  For  more  information  see,  Phys.  Rev.  E,  vol.  65,  pp.  056607 
(2002). 

5.  Discrete  Ginzburg-Landau  solitons 

The  complex  Ginzburg-Landau  (GL)  equation  is  known  to  play  a  ubiquitous  role  in 
science.  This  equation  is  encountered  in  several  diverse  branches  of  physics,  such  as,  for 
example,  in  superconductivity  and  superfluidity,  nonequilibrium  fluid  dynamics  and 
chemical  systems,  nonlinear  optics,  and  in  Bose-Einstein  condensates.  In  optics,  the 
discretized  version  of  the  Ginzburg-Landau  equation  is  often  used  to  model 
semiconductor  laser  arrays  or  arrays  of  highly  nonlinear  semiconductor  amplifiers 
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(SOAs).  We  have  demonstrated  (physical  Review  E,  vol.  67,  pp.  026606,  2003)  that 
discrete  solitons  are  possible  in  Ginzburg-Landau  lattices.  More  specifically,  we  have 
shown  that  this  system  exhibits  unique  complex  dispersion  properties  associated  with  a 
Brillouin  zone.  As  a  result,  the  discrete  dispersion  /diffraction  behavior  of  a  GL  lattice 
differs  substantially  from  that  encountered  in  conservative  arrays.  In  general,  two  new 
types  of  discrete  GL  solitons  can  exist  under  the  same  conditions.  These  solutions  are 
located  either  at  the  base  or  at  the  edge  of  the  Brillouin  zone  and  bifurcate  at  different 
values  of  the  linear  gain.  As  a  result  of  discreteness,  we  find  that  this  system  exhibits 
features  that  have  no  counterpart  whatsoever  in  either  the  continuous  GL  limit  or  in  other 
conservative  discrete  models  as  in  discrete  nonlinear  Schroedinger  chains.  These  include, 
for  example,  on-site  and  intra-site  bright  DGL  solitons  that  can  both  be  stable  as  well  as 
new  bifurcation  types  that  cannot  be  identified  in  the  continuous  case. 
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CREOL  Research 


The  original  concept  for  the  MURI  was  based  on  using  collisions  between  multi- 
component  solitons  to  perform  computing  and  other  signal  processing  functions.  The  ratio 
of  the  components  changes  in  a  deterministic  way  on  collision  thus  leading  to  well-defined 
initial  and  final  states  which  are  interpreted  as  a  mathematical  operation. 

To  investigate  this  process  this  group  investigated  experimentally  and  theoretically: 

(a)  The  properties  of  spatial  solitons  in  Kerr,  quadratic  and  amplifying  media; 

(b)  Limitations  on  solitons  due  to  anisotropy  and  modulational  instability; 

(c)  Soliton  collisions  utilizing  different  types  of  solitons  and  their  applications; 

(d)  Other  unique  related  phenomena  and  their  potential  applications; 

Only  research  highlights  are  summarized  here.  Details  are  in  prior  reports. 

1 .  Quadratic  Media 

Quadratic  spatial  solitons  (QSS)  are  multi-component  (vector)  solitons  consisting  of 
all  of  the  frequency  components  coupled  via  the  second  order  optical  nonlinearity  x(2)  near 
a  wavevector  matching  condition.  A  variety  of  wavevector  matching  techniques  including 
Type  I  (birefringent)  and  quasi-phase-matching  (QPM)  by  periodic  poling  of  the 
nonlinearity  were  used,  primarily  in  second  harmonic  generation  (SHG)  geometries. 


FW:  0.2  GW/cm2 


FW:  12  GW/cm2 


1 . 1  Quadratic  Solitons  in  Bulk  Media 

Two  different  media  were  utilized  for  QSS  generation,  both  of  which  were  non- 
critically  phase-matched  (NCPM).  A  key  feature  of  NCPM  is  that  the  refractive  index 
curves  for  the  fundamental  (FW)  and  harmonic  (SH)  were  tangent  to  each  other  along 
crystal  axes.  Thus  the  angular  bandwidth,  i.e.  the  range  of  incident  light  angles  for  soliton 
generation  was  very  large,  very  useful  for  applications.  The  measured  bandwidth  of  the 
threshold  for  soliton  generation  versus  incidence  angle  for  Type  I  NCPM  in  KNb03  was  24 
degrees.  For  QPM  KTP  is  was  about  half  this  value,  both  being  about  two  orders  of 
magnitude  larger  than  for  Type  II  SHG. 

The  QSS  consists  of  in-phase  fundamental  and  harmonic  components.  This  phase 
sensitivity  makes  the  generation  of  QSS’es  by  injecting  both  field  components  too 
cumbersome  to  be  practical.  It  was  found  that  injecting  just  the  fundamental  and  allowing 

the  appropriate  second  harmonic  to  be  generated 
over  small  distances  led  to  well-formed  solitons  as 
shown  in  Figure  1(b).  In  principle,  QSS’es  can  be 
generated  with  variable  fundamental  to  harmonic 
field  ratio,  depending  on  the  soliton  power  which 
would  be  ideal  for  exchange  of  the  different  field 
components  during  a  collision.  Unfortunately,  using 
fundamental  only  excitation  crystal  optical 
anisotropy  and/or  small  input  beam  ellipticity  led  to 
the  generation  of  multiple  solitons  specifically 
oriented  in  space  and  with  powers  less  than  twice 
threshold,  too  limited  to  be  useful  for  computing. 
This  general  limitation  was  found  in  both  NCPM 
samples  investigated. 


FW:  20  GW/cm1!  FW:  Simulation!  SH:  20  GW/cn 


Fig.  1  (a)  Diffracted  low  power 
fundamental  beam,  (b)  QSS  just  above 
threshold,  (c)  Excess  radiation  emitted  at 
twice  threshold,  (d)  Fundamental  and  (f) 
harmonic  emitted  at  three  times  threshold, 
(e)  Simulation  for  cases  in  (d)  and  (e). 
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1.2  QSS  in  Slab  Waveguides 


In  parallel,  QSSes  were  also 
investigated  in  slab  (ID)  waveguides  made 
by  Ti:indiffusion  using  Type  I  and  QPM 
phase-matching.  Again  it  proved  possible 
to  generate  high  quality  QSSes  by  injecting 
only  the  fundamental  beam,  but  again 
multi-soliton  generation  occurred  at  high 
input  powers. 

The  QPM  samples  were  used  to 
investigate  for  the  first  time  the  ratio  of 
fundamental  to  harmonic  power  in  a  ID 
QSS  versus  its  spatial  width  for  a  fixed 
value  of  wavevector  mismatch.  The  results 
are  shown  in  Figure  2  and  are  in  good 
agreement  with  theory. 


1.3  Discrete  Quadratic  Solitons  :  ID 

Arrays  of  channel  waveguides  fabricated  in  sufficiently  close  proximity  that  their  fields 
overlap  constitute  a  weakly  coupled  array  known  as  a  ID  discrete  system.  Light  at  the 
fundamental  frequency  “diffracts”  throughout  the  array,  i.e.  spreads  from  input  channels 
outwards  through  neighboring  channels,  due  to  coupling  from  one  channel  to  another.  The 
second  harmonic  light  is  too  strongly  confined  to  each  channel  to  couple  between  them  and 
appears  in  channels  solely  due  to  conversion  from  the  fundamental  present  in  them.  When  the 
media  are  nonlinear,  discrete  solitons  can  be  generated  in  which  self-trapping  occurs  over  a 
discrete  number  of  channels,  one,  two  etc. 


x  k^d 


The  linear  dispersion  relation  was  investigated.  As 
predicted,  regions  of  normal  diffraction,  anomalous 
diffraction  separated  by  zero  diffraction  were  found  and 
verified  experimentally.  Such  arrays  will  be  useful  for 
electro-optic  beam  deflectors. 

Fig.  3  Discrete  dispersion  relation  for  the  propagation 
wavevector  kz  versus  the  relative  phase  between  adjacent 
channels  (kxd). 


We  were  able  to  observe  for  the  first  time 
both  staggered  (tt  phase  shift  between  fields  in 
adjacent  channels)  and  unstaggered  (in  phase 
fields)  discrete  quadratic  solitons.  Shown  in  Fig.  4 
is  the  unstaggered  case  and  clearly  the  agreement 
between  experiment  and  theory  is  excellent.  This 
was  also  the  case  for  the  staggered  solitons. 
Highly  confined  (single  channel)  solitons  were 
also  observed. 


Fig.  4  Intensity  profiles  of  the  diffracted  beam  (open  circles  -  calculated;  dashed  line  -  measured)  and 
a  discrete  unstaggered  soliton  tsolid  circles  -  calculated:  solid  line  measured!. 
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1 .4  Instabilities 


It  was  predicted  theoretically  that  noise  on  a  wide  beam  in  quadratically  nonlinear 
media,  either  continuous  or  discrete,  gives  rise  to  filamentation.  This  limits  the  utility  of 
quadratic  solitons  for  computing  or  signal  processing.  This  effect  was  investigated  in  both 
ID  slab  waveguides  and  in  ID  discrete  arrays. 

Shown  in  Fig.  5  is  the  evolution  with  increasing  input  power  of  a  fundamental  beam 
injected  into  a  QPM  Ti:indiffused  LiNbCb  slab  waveguide,  (a)  shows  diffraction  at  low 
input  powers  and  collapse  into  a  soliton.  As  the  power  is  increased  further,  (b),  side  lobes 
appear  in  the  pattern  and,  at  the  highest  power,  multiple  peaks  occur  which  change  with 
each  laser  shot,  the  signature  of  noise-induced  filamentation. 


Fig.  5  FW  profiles  .  (a)  dashed 
line  -  74  pm  input  beam;  dark 
line  -  diffracted  beam;  heavy 
dark  line  -  ID  QSS.  (b)  3-peaked 
line  -  output  at  twice  soliton 
threshold  power;  other  lines  - 
irreproducible  filamentation  at  3 

times  threshold 


Similar  effects  occurred  in  bulk  samples  in  which  case  the  irreproducible  patterns  occur  in 
2D.  This  was  observed  in  the  NCPM  Type  I  KNb03  sample. 

Modulational  instability  was  also  observed  in  the  discrete  arrays  with  results  similar 
to  those  that  will  be  discussed  later  for  AlGaAs  discrete  arrays. 


1.5  Collisions 

The  large  angular  bandwidth  discussed  previously  in  1 . 1  for  quadratic  solitons  in  bulk 
non-critical  phase  matching  (NCPM)  geometries  was  ideal  for  the  investigation  of  collisions 
between  quadratic  solitons.  Collisions  between  two  identical  solitons  in  KNbC>3  were 
investigated  in  both  coplanar  and  non-co-planar  geometries.  In  the  first  the  solitons  were  co- 
planar,  launched  at  a  relative  angle  of  0.17°.  The  details  of  the  collapse  into  a  single  soliton  at 
zero  phase  angle,  the  repulsion  at  a  phase  angle  of  n,  and  the  preferential  energy  exchange  at 
intermediate  angles  are  all  clear  in  the  images  obtained  at  the  output  of  the  crystal.  To  within 


phase 

5%,  energy  is  conserved  in  the  collision,  as  is  the  composition  of  fundamental  versus 
harmonic.  In  the  second  geometry,  the  solitons  propagate  in  planes  10pm  apart  with  a 
“collision  angle”  of  0.9°  when  projected  into  a  common  plane.  A  third  soliton  was  observed 
after  collision  over  a  range  of  relative  phase  angles.  These  results  indicate  that  QSS 
collisions  are  unsuitable  for  applications  where  the  relative  phase  may  be  ill-defined. 
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1.5  Related  Phenomena  and  Applications 
1.5.1  Switching 

Three  experiments  on  switching  of  beams  were  performed.  The  first  relied  on  the 
multi-soliton  generation  phenomenon  discussed  previously.  In 
brief  summary,  a  small  amount  (1%  intensity)  in  harmonic 
seeding  caused  the  two  subsidiary  solitons  in  Figure  1  to 
disappear.  This  is  one  of  the  rare  cases  where  a  weak  beam 
can  control  strong  beams. 

The  second  case  is  a  consequence  of  the  fundamental 
and  harmonic  field  components  of  the  soliton  being  in-phase. 

When  the  poling  pattern  in  PPKTP  is  shifted  by  half  a  period 
as  shown  in  Figure  7,  this  introduces  an  effective  barrier 

potential  which 
can  reflect  a 
soliton.  When  the 
incidence  angle  is 
increased,  a 

transition  from 
reflection  to 

transmission 

through  the  interface  occurs.  A  similar  transition 

occurs  as  the  power  is  increased.  This  can  be  used 
for  switching. 

The  general  idea  for  the  third 
application  is  summarized  in  Figure  9.  A 
strong,  SH  control  beam  (blue)  was  coupled 
into  a  selected  channel  of  an  array.  An  input 
signal  beam  traversed  the  array  at  the  non- 
diffracting  angle  and  interacted  with  the 
control  beam  to  produce  a  deflected  beam. 

For  strong  control  beams,  the  input  and  output 
signal  experienced  gain  in  this  parametric 
process.  For  input  signals  not  exactly  half  the 
control  frequency,  the  output  was  frequency 
shifted.  Demonstrated  was  the  generation  of  a  frequency  shifted  output  signal  (“idler”  at 


Fig.  8  Output  FW  and  SH 
soliton  components  for 
different  incidence  angles 
on  the  interface. 


Fig.  7  Sample-light  geometry. 


coni'101 
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Fig.  9  Interaction  geometry. 


~1 554nm)  from  an  input  signal  at 
optical  communications. 


T543.5nm.  This  can  be  used  for  frequency  shifting  in 


1.5.2  Discrete  T albot  Effect 

The  Talbot  effect  in  continuous  media  is  well-known  since  the  mid  1800s.  A  periodic 
pattern  is  reconstituted  in  space  every  integer  multiple  of  the  well-known  Talbot  distance.  The 
discrete  Talbot  effect  was  first  investigated  theoretically.  It  was  found  that  recurring  patterns 
only  occur  for  the  excitation  of  every  1st,  2nd,  3rd,  4th  and  6th  channel.  In  continuous  media 
there  is  no  such  restriction.  In  discrete  media  it  is  the  different  form  (periodic)  of  the  dispersion 
relation  found  there  that  is  responsible.  This  prediction  was  confirmed  in  a  series  of  low  power 
experiments  in  PPLN  arrays  at  1550nm. 
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2.  ID  Kerr  Solitons  in  AlGaAs  Below  Half  the  Bandgap 

Spatial  solitons  based  on  the  Kerr  effect  exist  in  only  ID  (slab  waveguides)  in 
continuous  media.  They  were  observed  both  as  scalar  and  vector  (Manakov,  two 
orthogonally  polarized  components)  in  the  1990s  and  their  properties  are  well-known. 

2. 1  Discrete  Kerr  Solitons 

The  early  pioneering  work  was  extended  to 
highly  localized  (essentially  single  channel)  and 
vector  (two  polarization)  discrete  Kerr  solitons. 

The  collapse  of  the  discrete  fields  into  a 
strongly  localized  (blocker)  soliton  with  increasing 
input  power  into  a  single  channel  is  shown  in 
Figure  10.  Despite  the  pulsed  nature  of  the  input 
fields,  the  collapse  is  strongly  reminiscent  of  a  first 
order  phase  transition.  Similar  results  were  obtained 
for  the  successful  excitation  of  vector  discrete  Kerr 
solitons.  This  is  very  different  from  the  gradual 
collapse  observed  in  slab  waveguides. 


All _ 211 _ a _ 211 _ Ml. 


Fig.  10  Collapse  of  the  output  intensity 
from  the  AlGaAs  array  versus  input 
power  into  the  central  channel. 


2.2  Collisions  and  Applications 

The  coherent  collision  between  two  ID  discrete  Kerr  solitons  was  investigated  to 
determine  its  suitability  for  computing  applications.  At  low  input  powers,  each  beam 
diffracts  and  interference  effects  similar  to  those  of  a  phased  array  antenna  occur  with 
changes  in  relativephase.  The  output  is  spread  typically  over  4-5  channels.  As  the  input 
power  is  increased,  the  output  collapses  into  single  channel  solitons  and  varying  the  phase 
causes  the  output  to  jump  from  channel  to  channel.  This  is  a  digital  beam  scanner.  At  very 
high  powers,  the  coupling  between  channels  disappears. 
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Fig.  11  (b)  Output  from 
array  as  a  function  of  the 
relative  phase  between 
two  discrete  solitons 
showing  phase-scanning 
of  single  solitons.  (d) 
High  power  output. 
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A  number  of  experiments  relevant  to  routing  in  two  dimensional  arrays  of  channel 
waveguides  were  also  performed.  They  involved  a  “blocker”  (strongly  localized  single 
channel)  soliton  interacting  with  a  beam  traversing  the  array  at  the  zero  diffraction  angle, 
i.e.  the  control  beam  does  not  spread.  Three  different  cases  were  investigated:  (1)  co¬ 
polarized  mutually  coherent  beams;  (2)  orthogonally  polarized  mutually  coherent  beams; 
and  (3)  orthogonally  polarized  mutually  incoherent  beams.  The  geometry  and 
representative  results  for  cases  2  and  3  in  (b)  and  (c)  are  shown  in  Figure  12.  Only  for  case 
(c)  was  the  deflection  of  the  blocker  phase-independent  and  hence  useful  for  routing.  Note 
the  deflection  is  digital.  Finally,  two  successive  control  beams  were  demonstrated  to  move 
the  blocker  solitons  successively  by  one  channel  each.  This  deflection  is  the  basic  operation 
of  a  router  switch. 
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Fig.  12  (l.h.s.)  Interaction  geometry  showing  the  blocker  soliton  (red),  the  control  beam  (green)  and  the 
partially  reflected  control  beam  (blue),  (middle)  The  phase  dependence  of  the  TE  polarized  blocker 
deflection  for  case  2,  TM  power  =  870W  -  incident  blocker  on  the  left,  (c)  Deflection  of  the  blocker  as 
a  function  of  TM  (control)  power  for  the  incoherent  beams  case  3  which  is  phase  independent. 


Low  Power  (0.4  jjlW  avg.) 


High  Power  (17jnWavg.) 


2.3  Filamentation  (Modulational  Instability) 

The  same  mechanisms  that  lead  to  solitons,  diffraction  balancing  self-focusing,  also 
lead  to  the  noise-induced  break-up  of  a  wide  intense  beam.  This  process  was  investigated  in 
both  a  ID  slab  and  in  a  ID  array  of  coupled  channel  waveguides. 

The  more  interesting  case  is  the  discrete  array  since  over  Vi  the  range  of  possible 
propagation  directions(7i/2>|kxd|)  these  two  mechanisms  lead  to  filamentation,  and  over  the 
other  half  (tc> | kxd | >tt/2)  no  fdamentation  (or  bright  solitons)  should  occur.  This  was 
predicted  theoretically.  The  experimental  results  are  shown  below  in  Figure  13. 

As  the  low  power  input 
wavefront  is  tilted,  the  beam 
traverses  the  array  at  different 
angles.  The  striations  are  due  to 
imperfections  in  the 
waveguides.  At  high  powers 
these  striations  narrow  into 
filaments,  indicative  of 
modulational  instability  gain. 
Note  that,  as  predicted,  there  are 
regions  in  which  filamentation 
does  not  (and  cannot)  occur. 
The  gain  coefficient  was  also 
measured  by  seeding  the  input 
beam  with  different  periods. 


kxD 


Fig.  13  Variation  with  input  wavefront  tilt  (kxD)  of  the  beam  at 
the  output  facet  for  two  different  input  power  levels,  high  and 
low.  Note  that  the  input  wavefront  tilt  causes  the  beam  to  move 
across  the  array.  The  weak  filamentation  at  low  powers  is  caused 
by  array  imperfections.  At  high  powers  strong  filamentation 
occurs  for  (jr/2>|kxd|  and  no  filamentation  for  7t>|kxd|>7i/2. 


3.  Amplifying  Media:  Semiconductor  Optical  Amplifiers  (SOAs) 


3.1  Dissipative  Solitons 

“Zero  parameter”  dissipative  spatial  solitons  were  observed  in  AlGaAs  semiconductor 
optical  amplifiers  (SOAs).  Diffraction  and  nonlinearity  are  balanced,  and  simultaneously  gain 
balances  loss.  For  a  specific  gain  value,  there  is  only  one  soliton  possible  with  a  specific  peak 
power  and  width,  i.e.  there  is  no  power-width  trade-off. 
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The  zero  parameter  properties  are  verified  in  Figure  14.  First,  there  is  no  transmission 

at  the  output  if  the  input  is  below 
threshold,  i.e.  it  is  all  absorbed. 
Above  the  threshold,  the  output 
(soliton)  light  intensity  remains  a 
constant,  independent  of  the  input. 
The  striations  were  verified  to  be 
due  to  imperfections  in  the  device 
since  they  were  fixed  in  position 
when  the  beam  was  translated 
transversely  across  the  sample. 
With  fixed  light  input  power,  there 
is  also  a  threshold  gain  for  soliton 
generation  and  the  output  intensity 
increases  with  the  gain.  These 
properties  are  all  in  excellent 
agreement  with  theory.  Note  that 
these  solitons  have  10’s  mW  peak 
powers!! 


Ifl  60  no  16D  3  *  J.s  *3  VB 

Input  intensity  (mW)  Current  (A) 


Fig.  14  Left-hand-side:  The  measured  light  intensity  at  the  output 
versus  the  input  light  power,  (a)  experiment,  (b)  theory. 
Right-hand-side:  The  measured  output  light  intensity  (fixed  input 
intensity)  versus  the  input  current,  (a)  experiment,  (b)  theory. 


Also  investigated  was  whether  two  component, 
incoherently  coupled  solitons  exist.  Two  cw  light  sources, 
TS  at  943nm  and  LD  (948nm)  were  simultaneously  and 
collinearly  injected  into  a  current  pumped  InGaAs  single 
quantum  well  SOA  and  the  output  of  each  beam  was 
monitored  as  a  function  of  input  TS  power.  The  results  are 
shown  in  Figure  15.  Clearly  the  two  beams  interact, 
exchange  energy  and  form  a  coupled  pair  by  the  output. 
The  beams  are  mutually  trapped  since  neither  beam 
diffracts.  The  results  are  in  excellent  agreement  with 
simulations  based  on  the  equations  previously  used  to 
predict  dissipative  solitons. 

Detailed  numerical  simulations  show  that  for  very 
large  propagation  distances  the  power  is  all  converted  into 
one  wavelength.  For  stability  the  zero  parameter  property 
requires  that  both  the  gain  and  the  diffraction  must  be 
identical  at  the  two  wavelengths,  otherwise  one  beam  ends 
up  with  all  of  the  input  power.  Note  that  the  widths  of  the 
two  beams  maintain  their  sech-like  shapes  and  their  same 
width  (self-similar  beams)  on  propagation.  Hence  these 
coupled  beams  are  not  really  “solitons”.  However,  over  the 
sample  lengths  involved,  they  are  sufficiently  locked 
together  to  be  potentially  useful  for  signal  processing. 

It  has  been  demonstrated  that  optical  discrete  solitons  are  possible  in  dissipative  arrays 
or  lattices.  Such  dissipative  lattices  can  be  implemented  by  using  periodically  patterned 
semiconductor  optical  amplifiers  and  absorbers.  Unlike  their  continuous  counterparts,  this  new 


TS  power  (mW) 

Fig.  15  Output  of  (a)  the  TS  beam 
and  (b)  the  LD  beam  as  the  TS 
power  is  increased. 
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class  of  solitons  exhibits  several  interesting 
characteristics  because  of  discreteness.  For 
example,  these  new  solutions  no  longer  belong 
to  a  zero  parameter  family  and  their  nature 
dpends  on  their  Floquet-Block  band  structure. 
Experiments  are  underway. 


3.2  Instabilities 


Fig.  17  MI  gain  versus  noise 
periodicity  at  three  different  pump 
current  levels  n. 


The  theory  of  modulational  instability  (MI)  was 
developed  for  dissipative,  gain  media  under  conditions 
appropriate  to  regions  for  which  solitons  have  been 
observed  previously.  It  was  found  that  fluctuations  in  the 
phase  and/or  amplitude  of  the  optical  field  or  the  carrier 
density  can  lead  to  filamentation. 

The  variation  in  the  instability  gain  coefficient 
with  the  spatial  frequency  “p”  of  the  noise  resembles 
closely  that  reported  previously  for  other  soliton  systems, 
as  shown  in  Figure  17.  There  are,  however,  some 
interesting  differences.  For  example,  the  gain  saturates 
with  increasing  pumping  current,  a  feature  not  found  in 

any  other  MI  system  in  the  small  signal  limit.  These 
calculations  also  yielded  insight  into  regimes  in  which 
solitons  with  sub-  10pm  widths  can  be  excited. 


3.3  Collisions 
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Fig.  18  Measured  (lower)  and  calculated 
(upper)  patterns  at  the  output  facet  for  two 
parallel  solitons,  20pm  wide  for  beam 
separations  of  22  pms  (optical  fields  overlap). 
Experimental  results  are  a  collage  of 
exneriments  with  variable  relative  nhase. 


The  interactions  between  these  solitons 
were  measured  or  two  different  input  separations. 
The  results  show  that  the  interaction  can  be  either 
local  (fields  overlap)  or  non-local  (no  field 
overlap).  In  the  local  case,  Figure  18,  the  beam 
scanning  with  variable  relative  phase  between  the 
input  solitons  is  very  clean  and  the  output  always 
corresponded  to  a  single  soliton  of  the  same  power 
as  the  input  power  of  each  individual  soliton.  This 
is  a  direct  consequence  of  the  zero  parameter 
nature  of  this  type  of  soliton. 

For  the  non-local  interaction,  the  gain  and 
index  change  accompanying  each  soliton  has  an 
effective  width  of  -50 pms  leading  to  a  non-local 
interaction.  Shown  in  Figure  19  is  the  evolution 
of  both  the  calculated  index  and  gain  distribution 
due  to  the  interaction  and  the  optical  fields.  Note 
that  a  third  maximum  occurs  in  the  index  change 
which  leads  to  and  supports  the  generation  of  a 
third  soliton.  The  initial  separation  can  optimized 
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to  give  three,  equal  power  solitons  at  the  output. 
This  third  soliton  has  been  observed  in  a  soliton 
crossing  interaction  (not  shown  here). 


Fig.  19  Calculated  field  intensity  (left-hand- 
side)  and  evolution  of  refractive  index  change 
with  distance  (right-hand-side)  for  two  input,  in 
nhase.  solitons. 


4.  Related  Phenomena:  Switching  in  a  Dual  Core  Photonic  Crystal  Fiber  (PCF) 

Via  a  collaboration  with  OFS,  we  have  obtained  a  prototype 
directional  coupler  fabricated  in  a  PCF,  i.e.  two  parallel  channel 
(defect  core)  device.  A  picture  of  the  cross-section  of  the  fibers  is 
shown  in  Figure  20.  Note  the  small  asymmetry  (larger  hole)  on  the 
right  hand  side  which  introduces  birefringence  between  the  two 
cores  in  this  prototype  fiber  and  limits  the  transfer  efficiency  for  this 
particular  fiber.  This  fiber  was  cut  to  different  lengths  and  the 
periodic  power  transfer  distance  between  the  two  cores  was 
measured  to  be  ~4mm. 

Power-dependent  transfer  the  two  cores  is  demonstrated  in 
Figure  20.  However,  the  plateau  observed  at  high  intensities  is 
spurious  because,  as  shown  in  Figure  20,  the  fraction  of  energy  transmitted  through  both  cores 
at  1550  nm  decreases  rapidly  with  increasing  input.  This  energy  is  lost  to  other  wavelengths 
due  to  multiple  nonlinear  optical  phenomena, 
control  beam  powers. 


Fig.  20  Cross-section  of  a 
dual  “defect  core” 
photonic  crystal  fiber. 
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Fig.  20  Left-hand-side:  Outputs  from  the  crossover  PCF  core  when  the  input  core  is  excited  and  the 
input  power  is  increased  for  light  polarization  parallel  and  perpendicular  to  the  line  joining  the  cores. 
Right-hand-side:  The  total  throughput  from  both  cores  at  1550nm  as  a  function  of  input  intensity. 
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Research  Highlights  -  Segev  and  Christodoulides 


1.  Solitons  in  periodic  structures:  Lattice  (“Discrete”)  Solitons 

1.1  Observation  of  two-dimensional  lattice  solitons 

In  March  2003,  our  team  has  made  an  important  breakthrough,  in  demonstrating  the  first 
solitons  in  a  two-dimensional  nonlinear  photonic  lattice  [Fleischer  et  al,  Nature,  vol.  422,  pp. 
147,  2003;  see  also  commentary  about  this  work  in  the  April  issue  o/Physics  Today,  2003,  in 
the  Physics  Update  section ].  This  experiment  was  in  fact  the  first  soliton  ever  demonstrated  in 
any  2D  periodic  structure  in  nature.  We  first  demonstrated  bright  2D  lattice  solitons  in  their 
simplest  realization:  in-phase  solitons  at  the  base  of  the  first  Brillouin  zone.  Then,  we  proceeded 
to  demonstrate  the  first  2D  spatial  gap  solitons  (self-trapped  “staggered”  wavepackets  arising 
from  the  edge  of  the  first  Brillouin  zone  in  defocusing  systems).  The  formation  of  the  2D 
nonlinear  photonic  lattice  relied  on  an  optical  induction  technique  in  which  a  2D  array  of 
waveguides  is  induced  in  a  nonlinear  medium.  The  waveguide  array  is  induced,  in  real  time,  in 
a  photosensitive  material  by  interfering  two  or  more  plane  waves.  A  separate  ‘probe’  beam  is 
launched  into  the  periodic  waveguide  array,  where  it  exhibits  discrete  diffraction  and,  at  a 
sufficiently  high  nonlinearity,  forms  a  lattice  soliton.  For  this  system  to  work,  it  is  essential  that 
the  waveguides  are  as  uniform  as  possible,  implying  that  the  interference  pattern  (inducing  the 
lattice)  must  not  change  in  the  propagation  direction.  For  this  to  happen  in  the  nonlinear 
medium,  the  interfering  waves  themselves  should  not  be  affected  by  the  nonlinearity.  At  the 
same  time,  the  probe  (soliton-forming)  beam  must  experience  the  highest  possible  nonlinearity. 
A  photoreffactive  material  with  a  strong  electro-optic  anisotropy  allows  this  scenario;  the 
interfering  beams  are  polarized  in  a  non-electro-optic  direction  and  the  probe  is  polarized  along 
the  crystalline  c  axis.  In  this  arrangement,  the  interfering  beams  will  propagate  mostly  linearly, 
while  the  signal  beam  will  experience  both  a  periodic  potential  and  a  significant  (screening) 
nonlinearity.  Our  method  of  optical  lattice  induction  allows  for  dynamic,  reconfigurable  arrays 
of  almost  any  geometry.  Typical  experimental  results  are  shown  below. 

Since  this  paper  appeared  (March  2003),  it  has  been  cited  more  than  110  times,  being 
among  the  top  1%  most  cited  in  optics  in  the  past  5  years.  The  ideas  and  methods  presented  in 
that  paper  are  currently  being  followed  by  many  groups  world-wide. . _ 
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Experimental  observation  of  the  first  2D  lattice  soliton  [figure  taken  from  Fleischer  et  al, 
Nature,  vol.  422,  pp.  147,  2003],  The  figure  presents  the  propagation  of  a  probe  beam  launched 
into  a  single  waveguide  at  normal  incidence  with  respect  to  the  2D  photonic  lattice. 

Upper  left:  Intensity  structure  of  the  probe  beam  at  the  lattice  output,  displaying  diffraction- 
broadening  at  low  nonlinearity  (200V). 

Upper  right:  Intensity  structure  of  the  probe  beam  at  the  lattice  output,  displaying  an  on-axis 
lattice  soliton  at  high  nonlinearity  (800V). 

Lower  left:  Reduction  of  probe  intensity  by  x8,  at  the  same  voltage  (800V),  results  in  recovery 
of  the  diffraction-broadening  pattern. 

Lower  right:  Interferogram,  showing  constructive  interference  of  peak  and  lobes  between  the 
soliton  and  a  plane  wave. 

1.2  Spatial  gap  solitons  and  higher  bands  solitons 

In  the  3  years,  our  team  has  studied  various  manifestations  of  spatial  gap  solitons,  which  are 
the  spatial  analog  of  the  temporal  gap  solitons  forming  when  light  interacts  resonantly  with  a 
nonlinear  periodic  system.  In  2003,  we  have  made  the  first  experimental  observation  of  a  bright 
spatial  gap  soliton  [Fleischer  et  ah,  Physical  Review  Letters,  vol.  90,  pp.  23902,  2003].  It  was 
done  in  a  ID  configuration,  in  an  optically-induced  lattice  similar  to  the  2D  lattice  described  in 
Section  1.1,  and  two  months  later,  we  have  reported  the  observation  of  the  first  2D  spatial  gap 
solitons  [Fleischer  et  ah,  Nature,  vol.  422,  pp.  147,  2003], 


1.3  Vortex-ring  “discrete”  solitons  in  2D  photonic  lattices 

The  observation  of  solitons  in  a  two-dimensional  nonlinear  photonic  lattice  (described  in  the 
section  above)  has  opened  up  many  new  fascinating  possibilities.  One  of  them  is  the 
observation  of  vortex  solitons  in  photonic  lattices.  In  a  later  study  (Fleischer  et  al,  Physical 
Review  Letters  vol.  92,  pp.  123904,  2004),  our  team  has  presented  the  experimental 
observation  of  both  on-site  and  off-site  vortex-ring  lattice  solitons  of  unity  topological  charge  in 
a  square  lattice.  The  experiments  utilized  the  optical  induction  technique  in  a  photoreffactive 
crystal,  for  which  the  nonlinearity  is  inherently  saturable.  The  experiments  were  accompanied 
by  a  theoretical  study  of  such  vortex-ring  lattice  solitons:  we  found  the  wave  functions  of  on-site 
and  off-site  vortex  lattice  solitons,  and  demonstrated  their  stability.  These  vortex  solitons  are 
fundamental  wave  structures  on  2D  nonlinear  lattices  and  can  thus  provide  valuable  insight  into 
other  systems  in  nature  where  similar  dynamics  can  be  potentially  observed.  For  example, 
matter  waves  in  optical  traps  and  charge  density  waves  in  planes  of  conductive  polymers  can 
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support  solitons  with  topological  charge.  All-optical  experiments  also  hold  much  potential  for 
application,  e.g.  a  2D  array  of  fiber  lasers  can  excite  a  vortex  mode  whose  output  can  drive  a 
phase-sensitive  object.  Typical  experimental  results  are  shown  below.  A  picture  from  this 
work  also  featured  on  the  cover  of  the  relevant  issue  of  Physical  Review  Letters  [Volume 
92,  number  12,  March,  2004]. 


Experimental  observation  of  vortex  lattice  solitons  [figure  taken  from  Fleischer  et  al, 
Physical  Review  Letters  vol.  92,  pp.  123904,  2004],  The  photographs  show  results  taken  at  the 
output  of  the  photonic  lattice. 

(a,b)  Linear  diffraction  results:  (a)  intensity  and  (b)  phase  information  formed  by  interference  of 
output  with  a  plane  wave. 

(c)  Typical  on-site  lattice  soliton  intensity  and  (d)  soliton  phase  from  an  interferogram. 

(e)  Typical  off-site  lattice  soliton  intensity  and  (f)  soliton  phase  from  an  interferogram 


1.4  Second-hand  vortex-ring  solitons  in  2D  photonic  lattices 

The  observations  of  2D  lattice  solitons  (described  in  section  1.1)  and  of  second-hand 
lattice  solitons  (described  in  previous  sections),  have  naturally  raised  the  possibility  of  observing 
higher-band  solitons  in  2D  photonic  lattices.  Indeed,  our  team  has  predicted  [Manela  et  al., 
Optics  Letters,  vol.  29,  pp.  2049,  2004],  and  later  on  experimentally  observed  [Bartal  et  al., 
Physical  Review  Letters,  2005]  second-hand  solitons  in  2D  photonic  lattices.  These  solitons 
arise  from  the  X  symmetry  points  of  the  2nd  band,  and  they  reside  in  the  gap  between  the  1st  and 
the  2nd  bands  of  a  square  lattice,  i.e.  they  are  2D  spatial  gap  solitons  with  vorticity.  These 
bright  lattice  solitons  are  characterized  by  a  ring-shaped  intensity  profile,  and  a  unique  phase 
structure  resembling  a  counter-rotating  vortex  array.  Surprisingly,  even  though  these  solitons 
possess  a  complex  structure  (especially  their  vortex-array  phase),  they  can  be  excited  by 
launching  a  simple  highly-localized  vortex-ring  beam. 
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Our  simulations  and  experiments  reveal  that,  under  proper  nonlinear  conditions,  a  simple 
vortex  ring  excitation  naturally  evolves  into  the  2nd-band  vortex-ring  lattice  soliton,  acquiring  its 
unique  counter-rotating  vortex-array  phase  structure  during  propagation.  At  low-intensities, 
such  2nd-band  excitations  display  intriguing  features  of  preferential  linear  diffraction  along  the 
lattice  axes,  which  stand  in  sharp  contrast  to  the  diffractive  behavior  of  lst-band  vortex 
excitations  which  takes  on  the  (square)  symmetry  of  the  lattice.  We  demonstrated  these  2nd-band 
vortex-ring  lattice  solitons,  and  studied  their  evolution  dynamics  experimentally  and 
numerically.  The  figure  below  shows  typical  experimental  results  with  such  solitons.  We 
anticipate  that  the  dynamics  observed  here  will  appear  in  other  nonlinear  periodic  systems,  such 
as  nonlinear  fiber  bundles,  photonic  crystal  fibers,  and  Bose-Einstein  condensates  in  the  near 
future. 


Experimental  observation  of  second-band  vortex  lattice  soliton.  Intensity  (a), 
power  spectrum  (b)  and  phase  (c)  of  the  excitation  beam  at  the  input  plane  of  the 
lattice,  (d)  Linear  diffraction  of  a  low-intensity  beam,  (e)  Intensity  distribution  of  the 
second-band  vortex  soliton,  along  with  (f)  its  Fourier  power  spectrum  with  respect  to 
“comers”  of  the  1st  Brillouin  zone,  (g)  Phase  information  of  the  2nd  band  vortex  soliton, 
obtained  by  interference  with  plane  waves,  (h)  Numerical  validation  of  the  phase 
information.  This  figure  was  taken  from  Bartal  et  al.,  accepted  to  Physical  Review 
Letters,  2005. 


1.5  Brillouin-zone  spectroscopy  of  linear  and  nonlinear  photonic  lattices 

During  the  past  year,  our  team  has  developed  a  novel,  real-time,  experimental  technique  for 
linear  and  nonlinear  Brillouin  zone  spectroscopy  of  photonic  lattices  [see  Bartal  et  ah,  Physical 
Review  Letters,  vol.  94,  pp.  163902,  2005],  The  method  relies  on  excitation  with  random- 
phase  (partially- incoherent)  waves  and  far-field  visualization  of  the  spatial  spectrum  of  the  light 
exiting  the  lattice.  This  method  allows  the  characterization  of  the  underlying  lattice  structure, 
while  mapping  out  the  borders  of  the  extended  Brillouin  zones  and  marking  the  areas  of  normal 
and  anomalous  dispersion  within  them.  Specifically  for  photonic  lattices  with  defects  (e.g., 
photonic  crystal  fibers),  the  technique  enables  far-field  visualization  of  the  defect  mode  (guided 
mode)  overlay  on  the  extended  Brillouin  zones  structure  of  the  lattice.  The  technique  is  general 
and  can  be  used  for  photonic  crystal  fibers  as  well  as  for  periodic  structures  in  areas  beyond 
optics. 
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The  ideas  facilitating  this  spectroscopy  method  rely  on  linear  and  nonlinear  propagation  of 
random-phase  waves  in  a  photonic  lattice.  To  map  out  the  momentum  space  of  a  lattice,  the 
lattice  should  be  probed  by  a  broad  spectrum  of  its  eigen-modes  in  the  mapped  region  of  k- 
space.  Preferably,  one  would  like  to  excite  the  lattice  with  a  broad  and  uniform  spectrum  of 
Bloch  modes,  residing  in  several  different  bands.  At  the  same  time,  the  probe  beam  must  also  be 
broad  enough  in  real  space,  sampling  a  large  enough  number  of  unit  cells  so  that  the  lattice  is 
well  defined.  It  is  therefore  desirable  to  probe  the  lattice  with  a  beam  that  is  broad  both  in  real 
space  and  in  k-space  simultaneously.  These  requirements  are  reconciled  simultaneously  by  a 
random-phase  probe  beam,  which  facilitates  homogeneous  excitation  of  several  Brillouin  zones 
with  a  beam  that  occupies  numerous  lattice  sites.  Moreover,  the  time-averaged  intensity  of  such 
an  incoherent  probe  beam,  which  excites  many  Bloch  modes,  has  no  speckles,  because  all 
interference  effects  (among  Bloch  modes)  are  washed  out  by  the  stochastic  fluctuations  upon 
the  incoherent  beam.  Thus,  one  can  probe  the  photonic  lattices  with  a  partially-spatially- 
incoherent  (random-phase)  beam,  which  has  a  uniform  spatial  power  spectrum  extending  over 
several  Brillouin  zones,  and  is  broad  enough  to  cover  numerous  lattice  sites.  Typical  results 
depicting  the  extended  Brillouin  zone  map  of  two  photonic  lattices  (a  square  lattice  and  an 
hexagonal  lattice)  are  shown  below. 


Linear  mapping  of  the  edges  of  Brillouin  zones  of  square  and  hexagonal  lattices  [figure 
taken  from  Bartal  et  al.,  Physical  Review  Letters,  vol.  94,  pp.  163902,  2005],  (a)  Interference 
pattern  of  the  array  waves  inducing  the  photonic  lattice  with  the  square  symmetry,  (b)  Fourier 
spectrum  of  the  incoherent  probe  beam  (broad  circle)  and  lattice  forming  beams  (four  sharp 
peaks)  at  the  input  to  the  lattice,  (c)  Calculated  extended  Brillouin  zone  map  of  a  square  lattice, 
(d)  Experimental  Fourier  spectrum  of  the  probe  beam  at  the  lattice  output,  (e)-(h)  Same  as  (a-d) 
with  a  hexagonal  lattice. 

The  Brillouin  zone  spectroscopy  method  also  allows  direct  visualization  of  power 
spectrum  of  waves  propagating  in  lattices  with  defects.  The  figure  below  shows  an  example 
with  a  hexagonal  lattice  with  positive  and  negative  defects.  The  upper  row  depicts  the  near-field 
photograph  of  the  lattice  without  a  defect  (a),  and  with  a  positive  (b)  and  a  negative  (c)  defect. 
The  lower  row  (d-f)  shows  the  corresponding  (far-field)  power  spectrum  of  the  probe  beam  in 
each  case.  For  the  positive  defect,  the  guided  modes  (bound  states)  arise  from  total-intemal- 
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reflection  of  states  occupying  the  central  region  of  the  first  Brilluoin  zone.  Consequently,  the 
far-field  of  these  guided  modes  is  clearly  apparent  as  a  wide  spot  in  the  center  of  the  first 
Brilluoin  zone  of  (e).  On  the  other  hand,  for  a  negative  defect  (for  which  the  average  refractive 
index  is  lower  in  the  guiding  region),  waveguiding  arises  solely  from  Bragg  reflections,  and  not 
from  total-intemal-reflections.  Consequently,  the  central  region  of  the  Brilluoin-zone  picture  is 
empty  [not  populated;  the  central  “hole”  in  (f)],  and  the  guided  modes  are  the  states  whose 
momentum  arises  from  the  edge  of  the  first  Brillouin  zone  [bright  concentric  ring  in  (f))]. 


Defect  modes  in  a  hexagonal  lattice  [figure  taken  from  Bartal  et  ah,  Physical  Review  Letters, 
vol.  94,  pp.  163902,  2005],  Near  field  of  the  hexagonal  lattice  (a)  with  no  defect,  (b)  with  a 
positive  defect,  and  (c)  with  a  negative  defect.  Power  spectrum  of  the  incoherent  probe  beam 
exiting  the  hexagonal  lattice  (d)  with  no  defect,  (e)  with  a  positive  defect,  and  (f)  with  a 
negative  defect. 


Finally,  when  the  photonic  lattice  is  nonlinear,  the  underlying  self-focusing  (or  self- 
defocusing)  interaction  among  the  Bloch  modes  results  in  energy  exchange  between  these 
lattice  eigen-modes.  This  interaction  facilitates  a  method  for  distinguishing  between  the  regions 
of  normal  and  anomalous  diffraction  (dispersion)  of  the  underlying  lattice,  and  mapping  out  the 
boundaries  between  dispersion  of  opposite  signs.  In  these  experiments,  the  probe  beam  is 
propagating  nonlinearly  in  the  photonic  lattice,  and  it  drives  the  nonlinear  interaction  by 
inducing  a  broad  defect  in  the  lattice  structure,  which  in  turn  causes  energy  exchanges  between 
Bloch  modes.  Typical  experimental  results  are  shown  in  the  figure  below.  The  excitation  (far- 
field)  power  spectrum  of  the  incoherent  probe-beam  is  shown  in  (a).  When  the  lattice  exhibits  a 
nonlinearity  of  a  self-focusing  type,  the  wide  probe  beam  induces  a  wide  (and  shallow)  positive 
defect  in  the  lattice  (deeper  potential  /  increased  refractive  index).  The  presence  of  such  a  defect 
causes  Bloch  modes  with  anomalous  diffraction  to  transfer  power  to  Bloch  modes  experiencing 
normal  diffraction.  This  is  depicted  in  (b)  where  the  waves  residing  at  the  regions  of  normal 
diffraction  have  a  considerably  higher  intensity  than  the  low-intensity  regions  of  anomalous 
diffraction.  In  a  similar  fashion,  when  the  lattice  exhibits  a  nonlinearity  of  the  self-defocusing 
type,  energy  transfers  from  normal-diffraction  Bloch  modes  to  anomalous-diffraction  states, 
highlighting  the  higher  intensity  (anomalous  diffraction)  regions  in  (c),  adjacent  to  the  lower 
intensity  (normal  diffraction)  regions  in  the  same  far-field  power  spectrum. 
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Dispersion  mapping  of  a  square  latltice  [figure  taken  from  Bartal  et  al.,  Physical  Review 
Letters,  vol.  94,  pp.  163902,  2005  ].  (a)  Fourier  spectrum  of  the  probe  beam  (circle)  and  of  the 
lattice-forming  beams  (four  dots)  at  the  input.  Experimental  Fourier  spectrum  of  the  probe  beam 
at  the  output  face  of  the  square  lattice  under  (b)  self-focusing  nonlinearity  and  (c)  self- 
defocusing  nonlinearity 

To  conclude  this  section,  the  technique  for  linear  and  nonlinear  Brillouin  zone  spectroscopy  of 
photonic  lattices,  albeit  being  published  last  month,  has  already  generated  much  interest.  The 
method  is  general  and  can  be  used  for  any  periodic  optical  structure,  including  photonic  crystals 
and  photonic  crystal  fibers.  The  magazine  Physics  Today  has  published  a  commentary  about 
this  work  in  the  May  2005  issue,  in  the  Physics  Update  section. 


1.6  Random-phase  lattice  solitons 


Nonlinear  systems  with  inherent  periodicity  are  abundant  in  nature  Examples  may  be 
found  in  biology,  chemical  physics,  nonlinear  optics,  Josephson-junction  ladders,  Bose-Einstein 
condensates  (BEC),  etc.  Optics  especially  has  provided  important  advances  and  holds  great 
promise  for  applications  (e.g.  nonlinear  photonic  crystals).  The  main  feature  of  wave 
propagation  in  periodic  systems  is  the  interference  of  waves  reflected  from  the  lattice,  which 
drives  the  propagation  dynamics.  These  interference  effects  depend  crucially  on  the  coherence 
of  the  waves.  However,  in  nonlinear  periodic  lattices,  only  coherent  waves  have  been 
considered  thus  far.  Since  most  waves  encountered  in  nature  are  only  partially  coherent,  theories 
assuming  perfectly  coherent  waves  are  idealizations,  which  are  nevertheless  accurate  when  the 
characteristic  length  describing  the  coherence  of  waves  exceeds  by  far  the  characteristic 
dimension  of  the  system  (e.g.,  the  lattice  spacing).  However,  when  the  two  length-scales  are 
comparable,  the  interference  effects,  and  consequently  the  dynamics,  will  depend  on  the 
interplay  between  the  statistical  (coherence)  properties  of  the  waves  and  the  lattice  periodicity. 

In  paper  published  in  2004  [Buljan  et  al.,  Physical  Review  Letters,  vol.  92,  pp.  223901, 
2004],  our  team  has  present  the  first  theoretical  study  on  the  propagation  of  partially  coherent 
waves  in  nonlinear  periodic  lattices,  and  predicted  the  existence  of  random-phase  solitons  in 
nonlinear  periodic  lattices.  These  solitons  are  found  in  media  where  the  nonlinear  response  time 
is  much  longer  than  the  characteristic  time  of  random  phase  fluctuations.  A  random-phase 
lattice  soliton  forms  when  the  time-averaged  intensity  of  an  incoherent  wavepacket  induces  a 
defect  in  the  periodic  potential  which  has  multiple  bound  states;  the  wavepacket  "binds"  itself  to 
that  defect  by  randomly  populating  these  states  in  a  self-consistent  fashion  (on  the  time- 
average).  For  random-phase  lattice  solitons  to  exist,  their  intensity  profiles,  power  spectra,  as 
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FB  spectum  FT  spectrum 


well  as  their  statistical  (coherence)  properties  must  conform  to  the  lattice  periodicity.  The  main 
theoretical  findings  are  summarized  in  the  theory  figure  below. 


Theory:  properties  of  random-phase 
lattice  solitons  -  RPLS 

(figure  taken  from  Buljan  et  al.,  Physical 
Review  Letters,  vol.  92,  pp.  223901,  2004) 
(a)  Propagation  constants  of  the  modal 
constituents  of  RPLS.  (b,c)  The  power 
spectrum  on  an  RPLS  in  the  Fourier  and  in 
the  Floquet-Bloch  bases,  respectively,  (d) 
Intensity  profile  of  an  RPLS  (black)  and 
the  corresponding  correlation  distance 
across  it.  (e)  the  mutual  correlation 
function  of  the  RPLS  example,  (f) 
Characteristic  diffraction-broadening  when 
(d)  in  the  absence  of  nonlinearity. 


The  experiments  followed  soon  thereafter,  and  earlier  this  year  our  team  has  published 
the  first  observation  of  random-phase  solitons  in  photonic  lattices  [Cohen  et  al.,  Nature,  vol. 
433,  pp.  500,  2005].  A  typical  observation  of  a  random-phase  lattice  soliton  (RPLS)  is  depicted 
in  the  figures  below. 


Experimental  images  at  the  input  of  the  2D  photonic  lattice  [figure  taken  from  Cohen  et  al., 
Nature,  vol.  433,  pp.  500,  2005], 

(a)  Real-space  picture  of  induced  square  lattice  and  the  random-phase  probe  beam,  (b)  Power 
spectrum  of  the  incoherent  probe  beam  (wide  circle)  and  lattice  beams  (four  dots),  which  define 
the  comers  of  the  first  Brillouin  zone,  (c)  Same  as  (b)  but  with  coherent  probe  beam  that  has  the 
same  size  in  real  space  (as  in  (a)). 
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Experimental  observation  of  random-phase  lattice  solitons  [figure  taken  from  Cohen  et  al., 
Nature,  vol.  433,  pp.  500,  2005  ]. 

Shown  are  images  at  the  output  of  the  2D  photonic  lattice,  (a)  Linear  lattice  diffraction  of  the 
incoherent  beam  at  low  intensity,  (b-d)  Real  space  pictures  of  random-phase  lattice  solitons 
centered  (b)  on  a  site,  (c)  between  two  sites,  and  (d)  between  four  sites.  The  insets  show  3D 
perspectives  of  the  solitons.  Note  that  (a)-(d)  have  the  same  scale,  (e)  Power  spectrum  of 
random-phase  lattice  solitons  displaying  its  square  symmetry  and  multi-humped  structure  with 
the  peaks  located  in  the  normal  diffraction  regions  of  the  first  two  bands.  The  four  dots 
correspond  to  the  power  spectrum  of  the  array  beams. 


1.7  Modulation  instability  and  spontaneous  formation  of  solitonic  patterns  with 

incoherent  “white”  light 

One  of  the  most  fundamental  phenomena  in  nonlinear  dynamics  is  the  spontaneous 
formation  of  ordered  structures,  known  as  pattern  formation  or  modulational  instability  (MI) 
and  it  has  been  observed  and  studied  throughout  all  fields  of  science  such  as  fluids  dynamics 
and  plasmas,  biological  and  chemically  reacting  systems,  optics  and  many  others.  Essentially 
MI  is  a  process  in  which  a  nonlinear  system  in  a  spatially  homogeneous  state  is  becoming 
unstable  to  small  fluctuations  and  perturbations  with  a  specific  periodicity.  The  perturbations 
are  enhanced  through  the  nonlinearity,  until  diffraction  or  dispersion  counterbalances  the  effect 
and  a  regular  pattern  (typically  stripes,  rectangles  or  hexagons)  is  formed.  In  optics,  MI  has 
been  reported  in  many  materials  with  diverse  nonlinearities,  and  although  the  mechanisms 
driving  MI  are  different,  the  MI  process  itself  possesses  similar  properties  and  behavior.  Until 
recently,  pattern  formation  in  nonlinear  optical  systems  was  studied  with  fully  coherent  light, 
relying  on  the  intuition  that  spatial  correlations  are  necessary  for  wave  amplification.  However, 
5  years  ago,  we  have  shown  theoretically  and  experimentally  [Kip  et  al.,  Science,  vol.  290,  pp. 
495,  2000]  that  a  partially-coherent,  yet  quasi-monochromatic,  light  beam  can  also  undergo  MI 
breakup  resulting  in  a  periodic  array  of  spatially  incoherent  ID  or  2D  soliton-like  filaments. 
This  incoherent  MI  is  distinguished  from  coherent  MI  by  the  presence  of  a  distinct  threshold  for 
the  instability  to  occur.  That  is,  MI  occurs  only  when  the  nonlinear  self-focusing  is  strong 
enough  to  compensate  for  the  diffraction  of  light,  determined  by  the  spatial  correlation  function. 
Another  mandatory  condition  for  incoherent  MI  and  pattern  formation  to  occur  is  a  non- 
instantaneous  response  of  the  nonlinear  material,  so  that  the  nonlinear  change  of  the  refractive 
index  is  a  function  of  the  time-averaged  intensity  (and  not  of  the  instantaneous  intensity). 

Motivated  by  the  observation  of  optical  spatial  solitons  made  of  incoherent  white  light 
emitted  from  a  bulb  [Mitchell  and  Segev,  Nature,  vol.  387,  pp.  880,  1997],  our  team  has  started 
to  investigate  MI  of  white  light  (that  is,  light  which  is  both  spatially  and  temporally  incoherent). 
First,  a  new  general  theory  was  developed,  being  able  to  handle  the  propagation  of  incoherent 
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"white  light"  in  non-instantaneous  nonlinear  systems  [see  Buljan  et  al.,  PRE,  vol.  66,  pp.  35601, 
2002],  The  theory  predicted  that  MI  should  occur  with  white  light,  in  spite  of  the  temporal  and 
spatial  incoherence.  Similar  to  the  spatially-incoherent  yet  quasi-monochromatic  case,  the 
instability  occurs  only  when  the  nonlinearity  exceeds  a  threshold  value.  However,  with  white 
light  (wide  temporal  spectrum)  the  theory  has  indicated  that  the  MI  process  possesses  some  new 
features:  the  system  becomes  unstable  for  all  wavelengths  together,  at  the  same  nonlinearity 
value,  and  all  the  temporally  frequencies  that  make  up  the  temporal  spectrum  of  the  white  light 
collectively  participate  in  the  formation  of  the  pattern.  Therefore  the  temporal  spectrum  of  the 
evolving  perturbation  self-adjusts  as  the  beam  propagates  to  match  the  collective  MI 
phenomenon 

With  this  theoretical  backing  in  hand,  our  team  has  made  the  first  experimental 
observation  of  modulation  instability  and  spontaneous  pattern  formation  with  incoherent 
“white”  light  emitted  by  an  incandescent  light  bulb  [see  Schwartz  et  al.,  Physical  Review 
Letters,  vol.  93,  pp.  223901,  2004],  .  We  have  shown  that,  above  the  MI  threshold,  a  pattern 
consisting  of  white  light  stripes  forms,  and  that  this  pattern  exhibits  unique  dependence  on 
wavelength.  We  have  shown  experimentally  that  the  MI  threshold  value  is  the  same  for  all 
temporal  frequencies,  indicating  that  the  process  is  indeed  collective.  Increasing  the  nonlinearity 
beyond  the  threshold  results  in  the  enhancement  of  the  modulation  depth  of  the  periodic  pattern 
until  it  saturates.  The  figure  below  shows  typical  experimental  results.  This  work  has  also 
featured  in  the  December  2004  issue  of  Physics  Today  in  the  Physics  Update  section. 
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Intensity  distribution  at  the  exit  face  of  the  nonlinear  crystal  for  different  strengths  of 
nonlinearity,  showing  the  well  defined  threshold  for  white-light  MI.  [Figure  taken  from 
Schwartz  et  al.,  Physical  Review  Letters,  vol.  93,  pp.  223901,  2004). 

1.8  Theory  of  incoherent  white  light  solitons 

Our  team  has  theoretically  investigated  the  properties  of  solitons  made  of  temporally  and 
spatially  incoherent  light,  that  is,  white  light  solitons.  Their  evolution  dynamics  revealed  that 
the  spatio-temporal  coherence  properties  of  the  input  light  beam  change  in  a  characteristic 
fashion,  and  self-adjust  to  form  a  soliton.  We  identified  the  characteristic  features  of  the 
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temporal  power  spectrum  and  spatio-temporal  coherence  properties  of  white  light  solitons. 
More  specifically,  the  spatial  intensity  profile  of  light  within  their  bandwidth  is  wider  (less 
localized)  at  lower  frequencies,  and  narrower  at  higher  frequencies.  Furthermore,  the  spatial 
correlation  distance  (across  the  soliton)  is  always  larger  for  lower  frequencies,  and  shorter  for 
higher  frequencies.  We  have  studied  white  light  solitons  with  two  generic  types  of 
nonlinearities:  the  saturable  nonlinearity  and  the  Kerr  nonlinearity  when  both  have  a  non- 
instantaneous  response  (see  Buljan  et  al.,  Optics  letters,  vol.  28,  pp.  1239,  (2003)].  In  a  follow 
up  study  we  have  investigated  the  properties  of  white  light  incoherent  solitons,  analytically,  in  a 
logarithmically-saturable  nonlinear  media  (Buljan  et  al.,  PRE,  vol  68,  pp.  036607,  2003).  We 
found  closed-form  solution  representing  temporally  and  spatially  incoherent  solitons.  These 
incoherent  solitons  were  found  to  have  an  elliptic  Gaussian  intensity  profde,  and  elliptic 
Gaussian  spatial  correlation  statistics.  The  existence  curve  of  such  a  soliton  connects  the 
strength  of  the  self-focusing,  the  spatial  correlation  distance  at  a  particular  frequency,  and  the 
characteristic  width  of  the  soliton.  From  the  existence  curve  it  follows  that  these  solitons  exist 
only  when  the  spatial  correlation  distance  is  smaller  for  higher  frequency  constituents  of  the 
light. 


1.9  Grating-mediated- waveguides 

Optical  waveguides  are  widely  used  in  modem  optoelectronic  systems.  Optical  fibers  serve 
as  the  backbone  of  today’s  high-bandwidth  telecommunication  networks  while  many  system 
components  utilize  the  properties  of  waveguide  structures.  Thus  far,  three  physical  mechanisms 
for  optical  waveguiding  are  known:  Total  internal  reflection  (TIR),  "photonic  band  gap"  (PBG) 
waveguiding,  and  the  CROW  (Coupled  Resonator  Optical  Waveguide)  systems. 

Recently,  our  team  has  proposed  and  experimentally  demonstrated  a  new  method  for  optical 
waveguiding:  Grating-Mediated  Waveguiding  (GMW)  [Cohen  et  al.,  Physical  Review 
Letters,  vol.  93,  pp.  103902,  2004],  A  GMW  consists  of  a  structure  uniform  in  the  z-direction 
(propagation  direction)  in  which  waveguiding  is  driven  by  a  shallow  ID  grating  in  x,  which  is 
perpendicular  to  both  the  propagation  direction  z,  and  to  the  confinement  direction  y.  This 
grating  has  a  bell-shaped  or  a  trough-shaped  amplitude  in  y,  that  is,  in  the  direction  normal  to 
the  grating  wave-vector  (see  panels  a  and  b  in  the  figure  below).  Waveguiding  in  this  system 
occurs  when  two  angularly  symmetric  Bragg-matched  beams  are  incident  upon  the  grating,  with 
an  appropriate  relative  phase.  These  beams  are  simultaneously  Bragg-reflected  from  the  grating, 
and  are  jointly  guided  in  the  direction  normal  to  the  grating  wave-  vector.  In  the  case  where  the 
probe  beams  are  not  Bragg-matched,  or  if  their  relative  phase  is  improper,  no  waveguiding 
occurs  and  the  beams  diffract. 

We  proved  the  GMW  idea  theoretically,  and  demonstrated  its  concept  experimentally  in  a 
grating-mediated-waveguide  induced  in  a  photorefractive  SBN  crystal.  We  induce  a  trough- 
shape  GMW  structure  by  interfering  two  ordinarily-polarized  plane  waves  (inducing  an  index 
grating),  and  illuminating  the  central  interference  region  with  another  ordinarily-polarized 
beam,  which  is  very  narrow  in  y  and  is  made  incoherent  with  the  interfering  waves  [panel  c  in 
the  figure  below].  To  test  waveguiding,  we  launch  the  probe  beams  into  the  structure  and 
observe  their  propagation.  The  probe  beams  are  extraordinarily-polarized,  hence  they  “feel”  the 
waveguide  structure,  yet  they  are  made  sufficiently  weak  so  that  do  not  affect  the  induced  index 
structure.  Typical  experimental  results  are  shown  in  the  figure  below. 
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It  is  important  to  emphasize  the  difference  between  GMW  and  PBG  or  TIR  waveguides. 
In  PBG  waveguides,  the  reflections  are  in  the  plane  defined  by  the  propagation  and  confinement 
directions,  while  for  GMW,  the  grating  is  perpendicular  to  the  direction  of  confinement. 
Likewise,  in  TIR  structures,  the  light  is  confined  in  regions  of  higher  refractive  index  with  no 
Bragg  reflections,  whereas  for  GMW,  the  waveguiding  is  based  on  Bragg  reflections  and  all  the 
waveguide  “layers”  have  the  same  average  refractive  index.  This  last  argument  is  readily 
apparent  in  "trough"  GMW,  where  waveguiding  is  achieved  even  though  light  is  not 
concentrated  in  the  regions  of  the  peaks  of  the  refractive  index.  To  summarize,  we  have 
proposed  and  demonstrated  grating-mediated-waveguiding,  in  which  confinement  is 
perpendicular  to  the  grating  wavevector.  As  with  all  other  waveguiding  mechanisms,  GMW 
will  undoubtedly  find  use  in  a  variety  of  photonic  applications.  Moreover,  the  GMW  structures 
described  here  will  work  in  other  wave  systems  as  well,  such  as  acoustic  systems  and 
condensate  matter  waves  in  periodic  potentials. 
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Grating-mediated  waveguiding.  (a,b)  Schematic  drawings  of  GMW  with  a  bell  (a)  and  a  trough  (b) 
transverse  amplitude,  (c)  Experimental  scheme  for  inducing  a  trough-shape  GMW  in  photorefractives.  (d-h) 
Experimental  photographs  of  a  probe  beam  at  the  input  (d)  and  output  (e)  of  the  waveguide  providing  an 
experimental  proof  of  concept  for  GMW.  (f-h)  Diffraction  in  unguided  conditions:  a  homogeneous  medium 
(f),  when  the  beams  are  not  Bragg-matched  with  the  index  grating  (g),  when  the  beams  are  Bragg-matched 
with  the  grating  but  have  the  "wrong"  phase  relative  to  the  grating  (h).  The  intensity  in  each  panel  is 
normalized  to  the  peak  intensity.  This  figure  was  taken  from  Cohen  et  ah,  Physical  Review  Letters,  vol.  93, 
pp. 103902,  2004. 
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Scientific  Progress  and  Accomplishments 

This  report  covers  our  activities  and  lists  the  most  significant  accomplishments 
for  the  covered  period  2000-2005.  Our  activities  were  focused  on  developing  enabling 
technologies  to  create  important  building  blocks  for  an  all-optical  temporal-spatial  soliton 
gateway  and  its  peripherals.  Here  we  report  the  following  highlights: 

1)  All-optical  wavelength  conversion 

2)  Investigating  optical  spectral  bistability  for  all-optical  switching  and  optical  memory 

applications 

3)  Quantum  theory  of  Manakov  solitons 

4)  Multicomponent  gap  solitons  in  superposed  gratings 

5)  Quantum  phase  noise  reduction  in  soliton  collisions 

1)  All-optical  wavelength  conversion 

Solitonic  signal  processing  uses  the  fundamental  characteristics  of  optical  solitons 
to  achieve  multiple  functionalities.  We  made  a  significant  breakthrough,  both 
theoretically  and  experimentally,  in  building  a  gateway  with  all-optical  wavelength 
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conversion  capabilities  using  a  Sagnac  interferometric  loop  with  an  SOA  at  asymmetric 
position  .  We  filed  a  patent  disclosure  for  this  invention. 

All-optical  wavelength  conversion  is  one  of  the  most  important  functionalities  for 
optical  signal  processing,  which  has  found  wide  applications  for  ultra  short  optical  pulses 
/  temporal  solitons.  We  showed  that  the  structure  could  be  used  for  non-retum-to-zero 
signals  as  well. 


2)  All  optical  bistability  with  application  for  all-optical  memory 

In  November  2002,  during  the  review  process  we  had  a  live  demonstration  of 
optical  spectral  bistability  in  a  semiconductor  ring  laser,  in  which  hysteresis  of  the  laser 
spectrum  was  observed  in  response  to  an  optical  control  signal.  The  lack  of  a  viable 
technology  for  bit  level  random  access  optical  memory  is  a  limiting  factor  towards  the 
development  of  truly  intelligent  optical  networks.  Optical  bistability  in  lasers  with 
semiconductor  gain  media  is  a  promising  approach  toward  developing  bit  level  optical 
memories.  First  we  demonstrated  a  Bistable  Semiconductor  Fiber  Ring  Laser  (BSFRL) 
that  exhibits  tunable  bistable  wavelength  operation  near  the  1.55pm  region.  The 
corresponding  hysteresis  curves  are  shown  in  the  figure  below. 
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Fig.  The  two  bistable  states  of  the  BSFRL.  (a)  State  1  shows  a  narrow  spectral  peak  at 
1570.6  nm;  (b)  State  2  has  two  broad  peaks  at  1570.6  nm  and  1579.0  nm. 


3)  Quantum  theory  of  Manakov  solitons 

In  an  effort  to  determine  the  potential  of  communications  and  computation  systems  using 
Manakov  solitons,  the  Princeton  team  has  developed  a  quantum  theory  of  Manakov 
solitons.  The  theory  is  used  to  illustrate  the  vacuum-induced  fluctuations  on  Manakov 
soliton  propagation  and  collision.  This  approach  is  used  to  evaluate  the  performance  of 
systems  using  Manakov  solitons. 


4)  Multicomponent  gap  solitons  in  superposed  gratings 

We  investigated  the  stable  propagation  of  multicomponent  gap  solitons  (strictly  speaking, 
solitary  waves)  in  superposed  gratings,  non-uniform  structures  consisting  of  overlapped 
gratings.  Multicomponent  gap  solitons  are  comprised  of  multiple  pulses,  with  different 
carrier  frequencies,  corresponding  to  the  multiple  gaps  in  a  superposed  grating.  We 
analyzed  the  existence  and  stability  of  such  localized  waves  through  direct  numerical 
simulation  of  the  nonlinear  coupled-mode  equations  which  describe  this  physical  system. 
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Parts  (a-c)  of  the  figure  show  a  stationary  gap  soliton,  and  parts  (d-f)  demonstrate  a 
moving  gap  soliton.  Parts  (a,d)  are  for  one  component,  (b,e)  for  two  components,  and 
(c,f)  for  three  components.  For  the  two  and  three  component  solitons,  we  see  stable 
propagation  along  with  emission  of  dispersive  radiation.  Due  to  the  periodic  boundary 
conditions  used,  the  radiation  wraps  around  and,  in  the  case  of  parts  (b)  and  (c),  collides 
with  the  soliton.  This  perturbation  does  not  effect  the  propagation,  further  supporting  the 
stability  of  these  localized  waves. 


Fig.  Multicomponent  gap  soliton.  See  text  for  details. 


5)  Quantum  phase  noise  reduction  in  soliton  collisions 

We  discovered  a  new  physical  effect — that  phase  noise  is  reduced  by  a  soliton  collision, 
because  of  a  negative  correlation  between  phase  fluctuations  induced  by  self-phase 
modulation  (SPM)  and  cross-phase  modulation  (XPM).  We  then  showed  how  this  effect 
can  be  exploited  to  improve  quantum  nondemolition  (QND)  measurements.  The  effect 
itself  is  straightforward  to  implement  since  it  only  requires  operating  in  a  new  parameter 
regime,  and  hence  does  not  require  any  new  materials  or  detection  techniques.  We  show 
that  both  the  optimization  of  phase  noise  reduction  and  accurate  QND  measurements 
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favor  short  propagation  distances,  small  wavelength  separation  between  solitons,  and,  in 
general,  a  larger  probe  amplitude  with  respect  to  the  signal.  This  last  property  is 
particularly  surprising,  given  the  direct  correspondence  between  amplitude  and  SPM- 
induced  phase  noise,  a  characteristic  which  has  motivated  the  choice  of  smaller  probe 
than  signal  solitons  in  experiments.  Higher  fidelity  QND  measurements  may  lead  to  new 
applications  in  quantum  information  with  continuous  variables,  as  well  as  improve  QND- 
based  entanglement  and  various  applications,  including  entanglement  purification  and 
teleportation. 

The  overall  phase  fluctuations  are  measured  through  a  ratio  between  the  collision- 
induced  phase  noise  variance  with 

20 

respect  to  the  variance  in  the  absence 
of  collision.  This  ratio  is  plotted  in 
the  figure  as  a  function  of  probe  n 
amplitude,  where  the  signal  § 
amplitude  is  kept  fixed  at  1.  The  ~  15 
increase  in  the  variance  ratio  with  5 
increasing  propagation  distance,  .2 
shown  in  part  (a)  of  the  figure,  arises  §> 
due  to  SPM-induced  phase  noise,  §"  10 
which  tends  to  drown  out  the  phase  CL 
reduction  at  long  propagation 
distances.  In  addition,  the  variance 
reduction  favors  a  larger  probe.  The  ® 
ratio  can  also  be  minimized  for  small 
wavelength  separation  between 
probe  and  signal.  This  trend  is 
shown  in  part  (b)  of  the  figure,  in  <g- 
which  $z=5$,  and  also  demonstrates  o 
decreased  fluctuations  of  a  larger  ^ 
probe  amplitude  with  respect  to  the  ® 
signal. 
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We  predict  a  twofold  improvement  ^ 
in  the  error  variance  of  a  QND 
measurement  using  this  effect. 
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Verlag,  Berlin,  2001)  pp  37-60. 

2.  B.  Luther-Davies  and  G.I.  Stegeman,  Materials  for  Spatial  Solitons,  book  chapter  in 
Spatial  Solitons,  S.  Trillo  and  W.  Torruellas  editors  (Springer- Verlag,  Berlin,  2001)  pp  19- 
35. 

3.  W.  Torruellas,  Y.  Kivshar  and  G.I.  Stegeman,  Quadratic  Solitons,  book  chapter  in 
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Spatial  Solitons,  S.  Trillo  and  W.  Torruellas  editors  (Springer-Verlag,  Berlin,  2001)  pp  127- 
168. 

4.  G.I.  Stegeman,  Experiments  on  Quadratic  Solitons,  book  chapter  in  Proceedings  of 
NATO  Advanced  Research  Workshop  on  “Soliton  Driven  Photonics”,  A.D.  Boardman  and 
A.P.  Sukhorukov  editors,  (Kluwer  Academic  Publishers,  Holland,  2001),  pp  21-3. 

5.  L.  Friedrich,  R.  R.  Malendevich,  G.I.  Stegeman,  Instability  of  Fast  Kerr  Solitons  in 
AlGaAs  Waveguides  at  1.55  Microns,  book  chapter  in  Proceedings  of  NATO  Advanced 
Summer  Institute  on  “Soliton  Driven  Photonics”,  A.D.  Boardman  and  A.P.  Sukhorukov 
editors,  (Kluwer  Academic  Publishers,  Holland,  2001),  pp  317-320. 

6.  R.R.  Malendevich,  H.  Fang,  R.  Schiek,  G.I.  Stegeman,  Experiments  on  Seeded  and 
Noise  Initiated  Modulational  Instability  in  LiNbOs  Slab  Waveguides,  book  chapter  in 
Proceedings  of  NATO  Advanced  Summer  Institute  on  “Soliton  Driven  Photonics,  A.D. 
Boardman  and  A.P.  Sukhorukov  editors,  (Kluwer  Academic  Publishers,  Holland,  2001),  pp 
219-222. 

7.  T.  Carmon,  R.  Uzdin,  C.  Pigier,  Z.  H.  Musslimani,  M.  Segev  and  A.  Nepomnyashchy, 
Rotating  Propeller  Solitons,  Chapter  in  the  book  Proceedings  of  NATO  Advanced  Summer 
Institute  on  “Soliton  Driven  Photonics”,  A.D.  Boardman  and  A.P.  Sukhorukov  editors, 
(Kluwer  Academic  Publishers,  Holland,  2001) 

8.  M.  Segev  and  D.  N.  Christodoulides,  Incoherent  Solitons:  Self-Trapping  of  Weakly- 
Correlated  Wave-Packets,  Chapter  in  the  book  Optical  Spatial  Solitons,  Editors:  S.  Trillo 
amd  W.  Torruellas,  Publisher:  Springer-Verlag,  Berlin,  2001. 

9.  J.  B.  Smathers,  P.  Ballet,  Haeyeon  Yang,  C.  L.  Workman,  V.  R.  Yazdanpanah  and  G.J. 
Salamo,  Nucleation  Structure  and  Morphology  Instabilities  of  InAs  Self-Assembled  Islands, 
book  chapter  in  Quantum  Dots,  A.  D.  Andreev  editor,  in  the  series  Optoelectric  Properties  of 
Semiconductors  and  Superlattices,  in  press,  2002. 

10.  P.  R.  Prucnal,  V.  Baby,.  D.  Rand,  B.  C.  Wang,  and  I.  Glesk,  “ All-optical 
Processing  in  Switching  Networks",  IEEE  LEOS  Newsletter  16  (4)  13-14  (2002). 

1 1 .  G.I.  Stegeman,  “ All-Optical  Switching ”,  in  the  Handbook  of  Optics,  edited  by  M.  Bass, 
J.  M.  Enoch,  E.  W.  VanStryland  and  W.  L.  Wolfe  (McGraw  Hill,  New  York,  2001)  pp  21.1- 
21.9 

12.  V.  Baby,  L.  Xu,  B.  C.  Wang,  I.  Glesk,  and  P.  R.  Prucnal,  “ Optical  Spectral 
Bistability  in  a  Widely  Tuneable  Fiber  Ring  Laser  and  its  Applications  for  Bit-Lev 

el  Optical  Memory ,”  Resent  Research  Development  in  Electronics,  Vol.  1  -  2002,  S.  G. 
Pandalai,  Editor,  Transworld  Research  Network,  Trivandrum,  India,  pp.  123-131. 

13.  I.  Glesk,  D.  Rand,  and  P.  R.  Prucnal,  "TOAD:  Ultrafast  all-optical  device  for 
packet  routing  in  future  optical  networks ,”  Recent  Research  Development  in  Optics, 
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Vol.  3,  Research  Signpost,  Trivandrum,  India,  2003,  pp.  89-102,  ISBN:81-271-0028-5. 
Invited. 


14.  I.  Glesk,  B.  C.  Wang,  L.  Xu,  V.  Baby,  and  P.  R.  Prucnal,  “ Ultra-fast  all-optical 
switching  in  optical  networks ,”  Progress  in  Optics  Vol.  45,  E.  Wolf,  editor,  Elsevier 
Science,  Netherlands,  2003,  pp.  53-117.  Invited. 

15.  T.  Carmon,  H.  Buljan,  M.  Soljacic,  and  M.  Segev,  “ Pattern  formation  in  optical 
cavities  with  incoherent  light ”,  Optics  and  Photonics  News,  Special  Issue:  Optics  in 
2003,  vol.  14(12),  December  (2003). 

16.  Ladislav  Jankovic,  Hongki  Kim,  George  Stegeman,  Silvia  Carrasco,  Lluis  Tomer 
and  Mordechai  Katz  ,  " Self-reflection  and  routing  of  multicolor  solitons,"  Optics  & 

Photonics  News  (OPN)  "Optics  in  2003",  14(12):24,  (2003) 

17.  George  I.  Stegeman,  “Spatial  Beam  Instabilities  Due  to  Instantaneous  Nonlinear 
Mechanisms”,  book  chapter  in  proceedings  of  NATO  ASI/SUSSP56  titled  Ultrafast  Photonics, 
Ed.  Alan  Miller,  D.  T.  Reid  and  D.  M.  Finlayson,  (Institute  of  Physics  Publishing,  London, 
2004)  pp  103-21 

18.  G.  Assanto  and  G.  I .  Stegeman,  “ Nonlinear  Optics  Basics:  Cascading ”,  in  Encyclopedia 
of  Modem  Optics,  edited  by  Robert  D.  Guenther,  Duncan  G.  Steel  and  Leopold  Bayvel, 
(Elsevier  Press,  2004,  Oxford),  pp  207-12 

19.  E.  Ultanir,  G.I.  Stegeman,  D.  Michaelis,  C.H.  Lange,  and  F.  Lederer,  “Dissipative 
Solitons  in  Semiconductor  Optical  Amplifiers”,  Dissipative  Solitons:  Lecture  Notes  in  Physics, 
edited  by  Nail  Akhmediev  (Springer- Verlag,  Heidelberg,  2005),  661,  pp37-54 

20.  N.  K.  Eftemidis  and  D.  N.  Christodoulides,  book  chapter  “ Discrete  Ginzburg-Landau 
solitons”,  in  “Dissipative  Solitons ”  edited  by  Akhmediev,  Nail  and  Ankiewicz,  Adrian, 
Springer- Verlag,  2005. 

21.  T.  Coskun  ,  M.  Segev,  and  D.  N.  Christodoulides,  book  chapter,  “ Incoherent  Solitons ”, 
in  the  Encyclopedia  of  Nonlinear  Science,  edited  by  A.C.  Scott,  Routledge,  Dec.  2004. 

22..  E.  DelRe,  M.  Segev,  D.  Christodoulides,  B.  Crosignani,  and  G.  Salamo,  “ Photorefractive 
Solitons ”  ,  book  chapter  in  Photorefractive  Materials  and  Their  Applications  1,  edited  by  P. 
Gunter,  Springer  2005. 


6c  Papers  presented  at  meetings,  but  not  published  in  full  conference  proceedings 

281  Conference  papers,  seminars  etc  presented  relevant  to  MURI  activities 

Invited  (94) 
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1.  G.  I.  Stegeman,  “Prospects  for  Higher  Nonlinearities”,  (invited)  NSF  Workshop 
on  Networks  for  the  21st  Century,  Washington  DC,  December  2000 

2.  G.I.  Stegeman,  “Experiments  with  Quadratic  Solitons”,  (invited)  3  lectures  at 
NATO  ASI,  Swinoujscie,  Poland,  September  2000 

3.  G.  I.  Stegeman,  “Spatial  Solitons:  A  New  Frontier  in  Nonlinear  Optics”,  (invited) 

1  hour  tutorial,  QELS2000,  San  Francisco,  May  2000 

4.  I.  Glesk,  R.  Runser,  and  P.  R.  Prucnal,  "New  Trends  in  Optical  Communications", 
(invited)  Proceedings  of  Xllth  Czech-  Slovak-Polish  Optical  Conference  on  Wave  and 
Quantum  Aspects  of  Contemporary  Optics,  p.  45,  Velke  Losiny,  Czech  Republic, 
September,  2000. 

5.  G.  I.  Stegeman,  invited  talk,  Spatial  Solitons  and  Nonlinear  Beam  Dynamics  in 
Quadratically  Nonlinear  Media,  colloquium  at  the  University  of  Arkansas,  Fayetteville, 
October  2001. 

6.  R.  Schiek,  invited  talk,  Modulational  Instability  of  Plane-Wave  Eigenmodes  in 
Quadratic  Nonlinear  Media,  Progress  in  Electromagnetics  Research  Symposium,  Osaka 
Japan,  July  2001. 

7.  G.  I.  Stegeman,  invited  talk,  Modulational  Instability:  The  Role  of  Dimensionality  and 
Measurement  of  Instability  Gain  Coefficients  in  Quadratic  Media,  ICONO’2001,  Minsk, 
Belaruss,  June  2001. 

8.  G.  I.  Stegeman,  invited  talk,  Spatial  Modulational  Instability  in  Planar  Waveguides, 
NOMA  2001,  Cetraro,  Italy  May  2001. 

9.  G.  I.  Stegeman,  invited  talk,  Periodic  Spatial  Instabilities  of  High  Intensity  Beams  in 
Waveguides,  QELS  2001,  Baltimore,  May  2001. 

10.  G.  I.  Stegeman,  seminar,  Spatial  Modulational  Instability  in  Planar  Waveguides,  Dept. 
Physics,  Frederich  Schiller  Universitat,  Jena  Germany,  April  2001. 

11.  G.I.  Stegeman,  seminars,  Spatial  Solitons:  An  Experimental  Overview,  University  of 
Queensland,  Brisbane  Australia;  University  of  Melbourne,  Melbourne  Australia;  University 
of  New  South  Wales,  Sydney  Australia;  Australian  National  University,  Canberra,  Australia, 
February  2001;  University  of  Rome  III,  May  2001. 

12.  D.  N.  Christodoulides,  Incoherent  Spatial  Solitons,  invited  talk,  IMACS  Conference  on 
Nonlinear  Evolution  Equations  and  Wave  Phenomena,  Athens,  Georgia,  April  9-12,  2001. 

13.  D.  N.  Christodoulides,  Incoherent  Spatial  Solitons,  invited  talk,  paper  QThKl, 
CLEO/QELS  2001,  Baltimore,  Maryland,  May  6-11,  2001. 
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14.  D.  N.  Christodoulides  and  M.  Segev,  Incoherent  Spatial  Solitons,  invited  talk,  paper 
WN3,  IEEE  LEOS  Annual  Meeting,  La  Jolla  California,  November  11-15,  2001. 

15.  M.  Segev,  Incoherent  Solitons:  Self-Trapping  of  Weakly-Correlated  Wave-Packets, 
Invited  Talk,  International  Conference  on  Correlation  Optics  99,  Chemivtsy,  Ukraine,  May 
2001. 

16.  M.  Segev,  M.  Soljacic,  D.  Kip,  Z.  Chen,  S.  Sears,  E.  Eugenieva,  and  D.  N. 
Christodoulides,  Incoherent  Modulation  Instability:  Pattern  Formation  in  Weakly-Correlated 
Nonlinear  Wave  Systems,  paper  TuQ4,  invited  paper,  OSA  Annual  Meeting,  Long  Beach, 
California,  October  14-18,  2001. 


18.  M.  Segev  and  D.  N.  Christodoulides,  Incoherent  Spatial  Solitons,  invited  talk, 
CLEO/QELS  2002,  Long  Beach,  California,  May  19-24,  2002. 

19.  M.  Segev  and  D.  N.  Christodoulides,  Incoherent  Modulation  Instability:  Pattern 
Formation  in  Weakly-Correlated  Nonlinear  Wave  Systems,  Invited  Talk,  Optical  Society  of 
America  Annul  Meeting  (OSA  01),  Long  Beach,  California,  October  2001. 

20.  M.  Segev,  Incoherent  Solitons:  Self-Trapping  of  Weakly-Correlated  Wave-Packets, 
Invited  Talk,  International  Conference  on  Correlation  Optics  99,  Chemivtsy,  Ukraine,  May 
2001. 

21.  M.  Segev,  Progress  on  Optical  Spatial  Solitons  ,  Invited  Lecture,  Optical  Soliton 
Workshop  2001,  Orlando,  Florida,  March  2001. 

22.  M.  Segev,  Progress  on  Optical  Spatial  Solitons  ,  Invited  Lecture,  Optical  Soliton 
Workshop  2001,  Orlando,  Florida,  March  2001. 

23.  Lei  Xu,  Bing  C  Wang,  Varghese  Baby,  Darren  Rand,  Ivan  Glesk,  and  Paul  R.  Prucnal, 
"All-optical  data  format  conversion  between  RZ  and  NRZ  based  on  wavelength  converter," 
Lasers  and  Electro-Optics  Society,  2002.  LEOS  2002.  The  15th  Annual  Meeting  of  the  IEEE, 
Volume:  1 , 2002,  Page(s):  49  -50  paper  MF2  (invited). 

24.  P.  R.  Prucnal,  B.  C.  Wang,  and  I.  Glesk,  “Systems  Perspective  on  Optical 
Processing/Packet  Networks,”  paper  TuGl,  In  2002  Digest  of  the  LEOS  Summer  Topical 
Meeting:  All-optical  networking  Existing&  Emerging  Architecture  &  Applications,  Chateau 
Mont  Tremblant,  QC  Canada,  pp.  15-16,  15-17  July,  2002.  Invited  paper 

25.  G.  I.  Stegeman,  "Nonlinear  Beam  Dynamics  in  Media  with  Quadratic 
Nonlinearities",  OSA  Annual  Meeting,  Orlando,  Sept.  2002 
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26.  G.  I.  Stegeman  “Optical  Propagation  in  Quadratic  Nonlinear  Media:  Plane  Waves, 
Finite  Beams  and  Solitons”,  Nonlinear  Optics  Applications,  NOA'  2002,  Miedzyzdroje 
Poland,  Sept.  2002. 

27.  G.  I.  Stegeman,  “Beam  Dynamics  in  Nonlinear  Media”,  4  lectures,  Scottish 
Universities  Summer  Schools  in  Physics,  St.  Andrews  Scotland,  Sept.  2002 

28.  G.  I.  Stegeman,  “Overview  of  Spatial  Soliton  Experiments”,  Soliton  Workshop, 
Varenna,  Italy,  Aug.  2002 

29.  G.  I.  Stegeman,  “Spatial  Soliton  in  Periodically  Poled  KTP  (PPKTP)”,  Nonlinear 
Optics,  Maui  Hawai,  July  2002 

30.  G.  I.  Stegeman,  “Optical  Propagation  in  Quadratic  Nonlinear  Media:  Plane 
Waves,  Finite  Beams  and  Solitons”,  plenary,  Laser  Physics  Workshop  (LPHYS'02), 
Bratislava  (Slovak  Republic),  July  2002 

31.  M.  Segev,  D.  N.  Christodoulides,  M.  Soljacic,  Z.  Chen,  D.  Kip,  and  S.  Sears, 
“Pattern  formation  and  clustering  of  solitons  in  nonlinear  weakly-correlated  wave- 
systems”,  invited  talk  QThL3,  CLEO/QELS  2002,  Long  Beach,  California,  May  19-24, 
2002. 

32.  M.  Segev,  D.  N.  Christodoulides,  M.  Soljacic,  Z.  Chen,  D.  Kip,  S.  Sears,  and  T. 
Carmon,  “Pattern  formation  and  clustering  of  solitons  in  nonlinear  weakly-correlated 
wave-systems”,  invited  talk  ThD3,  NLO  02,  Maui,  Hawaii,  July  29- August  2,  2002. 

33.  G.I.  Stegeman,  invited  talk,  Nonlinear  Beam  Dynamics  in  Media  with  Quadratic 
Nonlinearities,  OSA  Annual  Meeting,  Orlando,  Sept.  2002. 

34.  G.I.  Stegeman,  invited  talk,  ID  QPM  Solitons  in  Homogeneous  and  Discrete  LiNbOs 
Waveguides,  Nonlinear  Optics  Applications,  NOA'  2002,  Miedzyzdroje  Poland,  Sept.  2002. 

35.  G.I.  Stegeman,  invited  talk,  Beam  Dynamics  in  Nonlinear  Media,  Scottish  Universities 
Summer  Schools  in  Physics,  St.  Andrews  Scotland,  Sept.  2002. 

36.  G.  I.  Stegeman,  invited  talk,  Spatial  Soliton  in  Periodically  Poled  KTP  (PPKTP), 
Nonlinear  Optics,  Maui  Hawai,  July  2002. 

37.  G.  I.  Stegeman,  invited  talk,  Spatial  Solitons,  Laser  Physics  Workshop  (LPHYS'02), 
Bratislava  (Slovak  Republic),  July  2002. 

38.  G.  I.  Stegeman,  “Real  Optical  Beam  Instability  and  Coherent  Spatial  Soliton 
Experiments  With  Non-Ideal  Samples  and  Non  Ideal  Sources”,  UCLA  Workshop  on 
Emerging  Applications  of  the  Nonlinear  Schrodinger  Equations,  Los  Angeles,  Feb. 

2003 
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39.  Joachim  Meier,  G.  I.  Stegeman,  S.  Aitchison,  R.  Morendotti,  and  Y.  Silberberg, 
“Nonlinear  Interactions  in  Discrete  Systems”,  (invited  paper)  Annual  LEOS  Meeting, 
Tucson,  Arizona,  October  2003,  paper  MC2 

40.  S.  Carrasco,  D.  V.  Petrov,  J.  P.  Torres,  L.  Tomer,  H.  Kim  and  G.  I.  Stegeman, 
“Tandem  Multicolor  Solitons”,  (invited  paper)  Annual  LEOS  Meeting,  Tucson, 
Arizona,  October  2003,  paper  WC1 

41.  S.  Aitchison,  R.  Morandotti,  D.  Mandelik,  H.  S.  Eisenberg,  Y.  Silberberg,  J. 
Meier,  and  G.  I.  Stegeman,  “Discrete  Spatial  Solitons  in  Kerr  Media”,  (invited  paper), 
Frontiers  in  Optics  (Annual  OS  A),  Tucson  Arizona,  October  2003,  paper  MZ2 

42.  W.  Sohler,  W.  Grundkutter,  J.  H.  Lee,  Y.  L.  Lee,  Y.  H.  Min,  V.  Quiring,  R. 
Ricken,  H.  Suche,  R.  Schiek,  T.  Pertsch,  F.  Lederer.  R.  Iwanow  and  G.  I .  Stegeman  , 
“All-optical  wavelength  conversion,  amplification  and  switching  in  periodically  poled 
Ti:LiNb03  waveguide  structures”,  (invited  paper),  Annual  LEOS  Meeting,  Tucson 
Arizona,  October  2003,  paper  TuSl 

43.  E.  Ultanir,  G.  I.  Stegeman,  Chr.  Lange  and  F.  Lederer,  "Dissipative  Propagating 
Spatial  Solitons  (Autosolitons)  in  Semiconductor  Amplifiers",  (invited  paper)  Topical 
Problems  of  Nonlinear  Wave  Physics  (NWP-2003),  Nizhny  Novgorod  Russia, 
September  2003 

44.  L.  Jankovic.  H.  Kim,  S.  Polyakov,  G.  I.  Stegeman,  S.  Carrasco,  L.  Tomer,  Chr. 
Bosshard,  P.  Gunter,  M.  Katz  and  David  Eger,  “Interactions  of  Quadratic  Spatial 
Solitons  in  Non-Critical-Phase-Matching  Geometries”,  (invited  paper),  Laser  Physics 
Workshop  (LPHYS'03),  Hamburg  Germany,  August  2003 

45.  G.  I.  Stegeman,  E.  Ultanir,  Chr.  Lange  and  F.  Lederer,  “Dissipative  Solitons  in 
Semiconductor  Optical  Amplifiers”,  (invited  paper)  Banfi  Memorial  Workshop,  Pavia 
Italy,  June  2003 

46.  George  I.  Stegeman,  Ladislav  Jankovic,  Hongki  Kim,  Sergey  Polyakov,  S. 
Carrasco,  L.  Tomer,  Chr.  Bosshard,  P.  Gunter,  M.  Katz  and  David  Eger,  “Spatial 
Solitons  in  Non-Critically-Phase-Matched  Crystals”,  (invited  paper)  NOMA’2003, 
Cetraro,  Italy  June  2003 

47.  Joachim  Meier,  George  Stegeman,  H.S.  Eisenberg,  Y.  Silberberg,  R.  Morandotti 
and  J.S.  Aitchison,  “Discrete  Beam  Interaction  in  Waveguide  Arrays”,  (invited  paper), 
IMACS’2003,  Athens  Georgia,  April  2003 

48.  J.  W.  Fleischer,  M.  Segev,  N.  K.  Eftemidis,  and  D.  N.  Christodoulides, 
“Observation  of  two-dimensional  discrete  solitons  in  optically-induced  nonlinear 
photonic  lattices”,  invited  paper,  paper  QThKl,  CLEO/QELS  2003,  Baltimore, 
Maryland,  June  1-6,  2003. 
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49.  D.  N.  Christodoulides,  N.  Efremidis,  J.  Fleischer,  and  M.  Segev,  “Discrete 
solitons”,  invited  talk,  NOMA,  Cetraro,  Italy,  June  8-13  (2003). 

50.  N.  K.  Efremidis,  J.  Hudock,  D.  N.  Christodoulides,  J.  W.  Fleischer,  and  M.  Segev, 
“Nonlinear  waves  in  periodic  lattices”,  invited  talk,  IMACS,  Athens,  Georgia,  April  7- 
10,2003. 

51.  M.  Segev,  J.  W.  Fleischer,  D.  N.  Christodoulides,  N.  K.  Efremidis,  and  T. 
Carmon,  “Discrete  solitons  in  optically-induced  lattices”  invited  talk,  European 
Conference  on  Lasers  and  Electro-optics,  Munich,  Germany,  June  22-27,  2003. 

52.  J.  W.  Fleischer,  M.  Segev,  N.  K.  Efremidis,  and  D.  N.  Christodoulides,  “Discrete 
solitons  in  optically- induced  photonic  lattices”,  invited  talk,  paper  MZ3,  OS  A  Annual 
Meeting,  Tucson,  Arizona,  October  5-9,  2003. 

53.  J.  W.  Fleischer,  O.  Cohen,  G.  Bartal,  T.  Schwartz,  H.  Buljan,  O.  Manela,  M. 

Segev,  N.K.  Efremidis,  J.  Hudock,  and  D.  N.  Christodoulides,  “Solitons  in  2D 
optically-induced  photonic  lattices”,  invited  talk,  paper  MCI,  IEEE/LEOS  2003, 
Tucson,  Arizona,  October  26-30,  2003. 

54.  G  J.  Salamo,  "Engineering  quantum  structures  and  their  behavior"  MRS,  Fall 
Meeting  2003(Invited) 

55.  I.  Glesk,  L.  Xu,  V.  Baby,  P.  Patel,  D.  Rand,  and  P.  R.  Prucnal,  “An  Optical 
Platform  for  Ultra-fast  Signal  Processing  and  Networks”  Europe-  U.S.  -Japan 
Symposium  on  Ultrafast  Photonic  Technology,  Makuhari  Messe,  Chiba,  Japan,  July  15, 

2003  (invited  paper). 

56.  G.  Bartal,  O.  Cohen,  H.  Buljan,  J.W.  Fleischer,  O.  Manela,  and  M.  Segev,  “Fourier-space 
imaging  in  nonlinear  photonic  lattices”,  invited  talk,  Frontiers  in  Optics,  Rochester,  New  York, 
October  2004. 

57.  G.  I.  Stegeman,  E.  Ultanir,  Chr.  Lange  and  F.  Lederer,  “Dissipative  Solitons  in 
SOAs  and  Their  Interactions”,  (invited  paper)  Photonics  and  Imaging  Initiatives 
Workshop,  Tucson  Arizona,  January  21-23  (2004) 

58.  G.  I.  Stegeman,  R.  Schiek,  D.  Christodoulides,  J.  Meier,  R.  Iwanow,  S.  Aitchison, 

R.  Morandotti,  Y.  Silberberg,  T.  Pertsch,  F.  Lederer,  Y.  Min  and  W.  Sohler,  “Discrete 
Solitons,  Their  Instabilities  and  Interactions  in  Quadratic  and  Cubically  Nonlinear 
Media”,  (invited  paper)  Winter  School  in  Nonlinear  and  Stochastic  Optics,  Tucson 
Arizona,  January  10-13  (2004) 

59.  G.I.  Stegeman,  J.  Meier,  S.  Aitchison,  R.  Morendotti,  Y.  Silberberg  and  G. 

Salamo,  “Spatial  Solitons”,  (colloquium)  Dept,  of  Physics,  University  of  Toronto, 

April  (2004) 
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60.  T.  Pertsch,  G.  I.  Stegeman,  R.  Schiek,  R.  Iwanow,  F.  Lederer,  Y.  Min  and  W. 
Sohler  “Transparent  Switching  in  PPLN  Waveguide  Arrays”,  (invited  paper), 
CLEO’2004,  San  Francisco,  May  2004 


61.  G.I.  Stegeman,  J.  Meier,  S.  Aitchison,  R.  Morendotti,  Y.  Silberberg  and  G. 

Salamo,  “Spatial  Solitons  and  Instabilities  in  Discrete  Kerr  Media”,  (invited  paper) 

NOA'  2004,  Warsaw  Poland,  June  2004. 

62.  Robert  Iwanow,  T.  Pertsch,  G.  I.  Stegeman,  R.  Schiek,  F.  Lederer,  Y.  Min  and  W. 

Sohler  “Discrete  Quadratic  Solitons  in  PPLN  Arrays”,  (invited  paper),  NOA'  2004, 

Warsaw  Poland,  June  2004. 

63.  G.  I.  Stegeman,  E.  Ultanir,  Chr.  Lange  and  F.  Lederer,  “Dissipative  Solitons  in 
Semiconductor  Optical  Amplifiers  and  Their  Interactions”,  (invited  paper)  ICTON 
2004,  Wroclaw  Poland,  July  2004 

64.  G.I.  Stegeman,  J.  Meier,  S.  Aitchison,  R.  Morendotti,  Y.  Silberberg  and  G. 

Salamo,  “Spatial  Solitons  in  Discrete  Nonlinear  Media”,  (invited  paper)  Laser  Physics 
Workshop  (LPHYS'04),  Trieste,  Italy,  July  2004 

65.  G.I.  Stegeman,  J.  Meier,  S.  Aitchison,  R.  Morendotti,  Y.  Silberberg,  G.  Salamo, 

R.  Iwanow,  T.  Pertsch,  R.  Schiek,  F.  Lederer,  Y.  Min  and  W.  Sohler  “Experiments  in 
Spatial  Solitons”,  (invited  paper)  Workshop  at  International  Centre  for  Mathematical 
Sciences  (ICMS),  Edinburgh  Scotland,  July  2004 

66.  D.  N.  Christodoulides,  “Discrete  solitons”,  invited  talk,  NLGW  workshop, 

Toronto,  Canada,  March  27,  2004. 

67.  M.  Segev,  J.  W.  Fleischer,  O.  Cohen,  H.  Buljan,  G.  Bartal,  T.  Schwartz,  O. 

Manela,  N.  K.  Efremidis,  J.  Hudock,  and  D.  N.  Christodoulides,  “Wave  propagation 
and  solitons  in  2D  photonic  lattices”,  invited  talk,  CLEO/IQEC  2004,  San  Francisco, 

California,  May  16-21,  2004,  paper  JTuD3. 

68.  D.  N.  Christodoulides,  “Discrete  solitons”,  invited  talk,  ICMS  Workshop  on 
Mathematical  Issues  in  Nonlinear  Optics,  Edinburgh,  UK,  July  18-24,  2004. 

69.  D.  N.  Christodoulides,  “Discrete  optical  solitons”,  invited  talk,  Nonlinear  Optics, 
Waikoloa-Hawaii,  August  2-6,  2004. 

70.  D.  N.  Christodoulides,  “Two-dimensional  discrete  solitons  in  optically  induced  nonlinear 
photonic  lattices”,  invited  talk,  CLEO/IQEC  2004,  San  Francisco,  California,  May  16-21, 
2004. 

71.  G.I.  Stegeman,  E.  Ultanir,  Chr.  Lange  and  F.  Lederer,  “Dissipative  Solitons  in 
SOAs  and  Their  Interactions”,  (invited  paper)  Photonics  and  Imaging  Initiatives 
Workshop,  Tucson  Arizona,  January  21-23  (2004) 
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72.  G.  I.  Stegeman,  Roland  Schiek,  Demetri  Christodoulides,  Joachim  Meier,  Robert 
Iwanow,  Stewart  Aitchison,  Roberto  Morendotti,  Yaron  Silberberg,  Thomas  Pertsch, 

Falk  Lederer,  Yoohong  Min  and  Wolfgang  Sohler,  “Discrete  Solitons,  Their 
Instabilities  and  Interactions  in  Quadratic  and  Cubically  Nonlinear  Media”,  (invited 
paper)  Winter  School  in  Nonlinear  and  Stochastic  Optics,  Tucson  Arizona,  January  10- 
13  (2004) 

73.  G.  I.  Stegeman,  E.  A.  Ultanir,  D.  Christoulides,  C.  H.  Lange  and  F.  Lederer,  “Dissipative 
Solitons  in  SOAs  and  Their  Interactions”,  (invited),  LEOS  Annual  Meeting,  Puerto  Rico, 
November  (2004) 

74.  G.  I.  Stegeman,  “What  are  Optical  Solitons  and  Why  Are  They  Interesting”,  School  on 
Optics  and  Photonics  (1  hr  student  tutorial)  RIAO/OPTILAS  2004,  Porlamar  Venezuela, 
October  2004 

75.  G.  I.  Stegeman,  J.  Meier,  S.  Aitchison,  Y.  Silberberg,  R.  Morandotti,  H.  Yang,  and  G. 
Salamo,  “Discrete  Spatial  Solitons  in  ID,  Kerr,  Optical  Arrays”,  (invited)  RIAO/OPTILAS 
2004,  Porlamar  Venezuela,  October  2004 

76.  R.  Iwanow,  R.  Schiek,  G.  I.  Stegeman,  T.  Pertsch,  F.  Lederer,  Y.  Min  and  W.  Sohler, 
"Discrete  Quadratic  Solitons  in  ID  Arrays",  (invited),  SIAM  Meeting  on  Nonlinear  Waves  and 
Coherent  Structures,  Orlando,  October  (2004) 

77.  W.  Sohler,  W.  Grundkotter,  J.  H.  Lee,  Y.  H.  Min,  V.  Quiring,  H.  Suche,  R.  Schiek, 

T.  Pertsch,  F.  Lederer,  R.  Iwanow,  and  G.  I.  Stegeman,  “All-Optical  Signal  Processing  in 
Periodically  Poled  LiNb03  Waveguide  Structures”  (invited,  given  by  Wolfgang  Sohler), 
ECOC'04,  Stockholm  Sweden,  September  (2004) 

78.  D.  N.  Christodoulides,  “Discrete  and  lattice  optical  solitons”  invited  paper,  SIAM 
meeting,  Orlando,  Florida,  October  2-5,  2004. 

79.  G.  I.  Stegeman,  “Dawn  of  Discreteness  in  Optics”,  Stoicheff  70’th  Birthday  Celebration 
(invited),  Toronto,  December  2005 

80.  G.  I.  Stegeman,  “Discrete  Linear  and  Nonlinear  Optics”,  Colloquium  Lehigh  University, 
October  2005 

8 1 .  George  I.  Stegeman,  Joachim  Meier,  Demetri  Christodoulides,  Roberto  Morondotti, 

H.  Yang,  G.  Salamo  and  M.  Sorel,  “Discrete  "Blocker"  Solitons  for  Routing  in  AlGaAs 
Arrays”,  Frontiers  in  Optics  (invited),  Tucson,  October  2005 

82.  G.  I.  Stegeman,  E.  Ultanir,  D.  Christodoulides,  F.  Lederer,  and  C.  H.  Lange,  “Dissipative 
Spatial  Solitons  in  Semiconductor  Optical  Amplifiers”,  (invited  paper)  NLGW’2005,  Dresden 
Germany,  September  2005 

83.  G.  I.  Stegeman,  “Linear  and  Nonlinear  Discrete  Optics”,  (plenary  talk),  SPIE 
International  Congress  on  Optics  and  Optoelectronics,  Warsaw  Poland,  August  2005 
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84.  G.  I.  Stegeman,  “Nonlinear  Beam  Propagation  in  Quadratically  Nonlinear  Materials”, 
(tutorial)  CLEO  2005,  Baltimore,  May  (2005) 


85.  G.  I.  Stegeman,  “Discrete  Spatial  Solitons”,  (tutorial),  ICONO  2005,  St.  Petersburg 
Russia,  May  2005 

86.  G.  I.  Stegeman,  E.  A.  Ultanir,  D.  Christoulides,  C.  H.  Lange  and  F.  Lederer,  “Dissipative" 
Spatial  Solitons  and  Their  Interactions  in  SOAs”,  (invited)  ECIO  2005,  Grenoble  France,  April 
2005 

87.  D.  N.  Christodoulides,  “All-optical  routing  using  discrete  solitons”,  invited  paper,  paper 
5735-07,  Photonics  West,  San  Hose,  California,  January  22-27,  2005. 

88.  D.  N.  Christodoulides,  “Discrete  solitons  in  waveguide  arrays  and  optical  lattices”, 
Keynote  Speech,  The  fourth  IMACS  international  conference  on  Nonlinear  Evolution 
Equations  and  Wave  Phenomena,  Athens,  Georgia,  April  11-14,  2005. 

89.  D.  N.  Christodoulides,  “Discrete  solitons  in  waveguide  arrays  and  optical  lattices”, 
invited  talk,  ISE1,  ICONO/LAT  2005,  St.  Petersburg,  Russia,  May  1 1-15,  2005. 

90.  B.  Freedman,  G.  Bartal,  M.  Segev,  D.  Christodoulides,  J.  Fleischer,  “Optical  quasi- 
crystals-properties  and  dynamics”,  invited  talk,  paper  QThKl,  CLEO/QELS,  Baltimore, 
Maryland,  May  22-27,  2005. 

91.  M.  Segev,  “Solitons  in  photonic  lattices”,  plenary  talk,  Conference  on  Nonlinear  Waves, 
Integrable  Systems,  and  their  Applications,  Boulder,  Colorado,  June  2005. 

92.  M.  Segev,  “Solitons  and  nonlinear  waves  in  photonic  lattices”,  plenary  talk  of  the 
Division  of  Quantum  Optics,  Annual  Meeting  of  the  German  Physical  Society  (DPG),  Berlin, 
Germany,  March  2005. 

93.  M.  Segev,  “Solitons  and  nonlinear  waves  in  photonic  lattices”,  invited  talk,  French-Israeli 
Symposium  on  Nonlinear  and  quantum  Optics  (FRISNO  2005),  Ein  Bokek,  Israel,,  February 
2005. 

94.  M.  Soljacic  E.  Lidorikis,  J.D.  Joannopoulos,  L.V.  Hau,  M.  Segev,  K.  Steiglitz,  and  C. 
Anastassiou,  “New  windows  of  opportunity  for  all-optical  information  processing:  spatial 
solitons,  EIT  in  micro-cavities,  etc.”,  invited  talk,  SPIE  Conference  on  Optoelectronics  2005, 
San  Jose,  California,  January  2005. 


Contributed  (187) 

1.  P.  V.  Studenkov,  M.  R.  Gokhale,  J.  Wei,  W.  C.  Lin,  I.  Glesk,  P.  R.  Prucnal,  and  S. 
R.  Forrest,  "Integrated,  Asymmetric  twin-waveguide  ultrafast  all-optical  switch,"  IEEE 
LEOS'OO,  paper  WCC4,  pp.  645-646,  Rio  Grande,  Puerto  Rico,  November  2000. 
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2.  R.  J.  Runser,  C.  Coldwell,  P.  Toliver,  I.  Glesk,  and  P.  R.  Prucnal,  "A  Practical 
All-Optical  Sampling  Technique  for  High  Bandwidth,  Low  Energy  Optical 
Communication  Signals,"  IEEE  LEOS'OO,  paper  THJ3,  pp.  770-771,  Rio  Grande, 
Puerto  Rico,  November  2000. 

3.  B.  Wang,  I.  Glesk,  Robert  J.  Runser  and  P.  R.  Prucnal,  "Tunable  Parallel  Optical 
Delay  Line,"  IEEE  LEOS'OO,  paper  WZ3,  pp.  621-622,  Rio  Grande,  Puerto  Rico, 
November  2000. 

4.  D.  Zhou,  B.  Wang,  I.  Glesk  and  P.  R.  Prucnal,  "Perfectly  Synchronized  Bit- 
Parallel  WDM  Data  Transmission  Over  a  Single  Optical  Fiber,"  The  2000  IEEE/LEOS 
Summer  Topical  Meetings,  Aventura,  Florida,  July  2000. 

5.  R.  J.  Runser,  I.  Glesk,  and  P.  R.  Prucnal,  "Experimental  Demonstration  of  a  1.5  ps 
Demultiplexing  Window  for  High  Speed  Optical  Networks  Using  a  Forward-Pumped 
Mach  Zehnder  TOAD",  CISS2000,  Proceedings  of  34th  Annual  Conference  on 
Information  Sciences  and  Systems,  Princeton  University,  NJ,  March  ,  2000 

6.  D.  Zhou,  I.  Glesk  and  P.  R.  Prucnal,  "A  novel  optical  pulse  width  management 
device,"  in  CLEO  Technical  Digest  (Opt.  Soc.  Am.,  Washington  DC,  2000),  pp.  133. 

7.  P.  Toliver,  R.  J.  Runser,  I.  Glesk,  and  P.  R.  Prucnal,  "Comparison  of  Three 
Nonlinear  Optical  Switch  Geometries,"  in  CLEO  Technical  Digest  (OSA,  Washington 
DC,  2000),  pp.  254-255. 

8.  H.  Fang,  R.  Malendevich,  R.  Schiek  and  G.  I.  Stegeman,  “Spatial  Modulational 
Instability  in  One-Dimensional  LiNbOs  Slab  Waveguides”,  QELS  2000  Digest,  San 
Francisco  CA,  May  2000,  paper  QThP2 

9.  L.  Friedrich,  R.R.  Malendevich,  G.I.  Stegeman,  J.M.  Soto-Crespo,  N.N. 
Akhmediev  and  J.S.  Aitchison,  “Are  spatial  solitons  of  both  polarizations  stable  in  Ken- 
slab  waveguides?”,  QELS  2000  Digest,  San  Francisco  CA,  May  2000,  paper  QThP5 

10.  R.  Schiek,  C.G.  Trevino-Palacios,  J.  Meier,  H.  Fang,  L.  Friedrich,  G.I.  Stegeman, 
K.R.  Parameswaran,  M.-H.  Chou  and  M.M.  Fejer,  “High-power  second-harmonic 
generation  experiments  in  periodically  poled  lithium  niobate  waveguides”  CLEO  2000 
Digest,  San  Francisco  CA,  May  2000,  CWH2 

11.  H.  Fang,  R.  Malendevich,  G.I.  Stegeman  and  R.  Schiek,  “Characterization  of 
spatial  modulational  instability  of  second-order  nonlinear  eigenmodes  in  a  slab 
waveguide”,  Topical  Meeting  on  “Nonlinear  Optics:  Materials,  Fundamentals  and 
Applications”,  Kaua’i-Lihue  Hawaii,  August  2000,  TuBl  1 

12.  R.R.  Malendevich,  H.  Fang,  R.  Schiek,  G.I.  Stegeman  “Experiments  on  Seeded 
and  Noise  Initiated  Modulational  Instability  in  LINBO3  Slab  Waveguides”,  NATO  ASI 
on  Soliton  Driven  Photonics”,  Swinoujscie,  Poland,  September  2000 
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13.  R.R.  Malendevich,  L.  Friedrich,  R.R.  Malendevich,  G.I.  Stegeman,  J.M.  Soto- 
Crespo,  N.N.  Akhmediev  and  J.S.  Aitchison,  “Instability  of  Fast  Kerr  Solitons  in 
AlGaAs  Waveguides  at  1.55  Microns”,  NATO  ASI  on  Soliton  Driven  Photonics”, 
Swinoujscie,  Poland,  September  2000 

14.  Yan-qing  Lu,  Min  Xiao  and  Gregory  J.  Salamo,  Wide-Bandwidth  High-Frequency 
Electro-optic  Modulator  Based  on  Periodically  Poled  LiNbOs,  CLEO  2001 . 

15.  Xiaodong  Mu,  Yujie  J.  Ding,  Haeyeon  Yang,  and  Gregory  J.  Salamo,  Cavity-Enhanced 
and  Quasi-Phase-Matched  Multi-Order  Reflection-Second-Harmonic  Generation  From 
GaAs/AlAs  and  GaAs/AlGaAs  Multilayers,  CLEO  2001. 

16.  R.  Uzdin,  M.  Segev,  G.  Salamo,  Self-Focusing  in  InP,  CLEO  2001 

17.  P.  Ballet,  H.  Yang,  C.  L.  Workman  and  G.J.  Salamo,  GaAs  Dots  on  GaP,  MRS,  2001. 

18.  Haeyeon  Yang,  P.  Ballet,  and  G.J.  Salamo,  Formation  of  Quantum  Wires  and  Dots  on 
InP (001)  by  As/P  Exchange,  MRS,  2001  . 

19.  Bing  C.  Wang,  Varghese  Baby,  Lei  Xu,  Ivan  Glesk,  Paul  R.  Prucnal,  A  Fast  Optical 
Switch  Based  on  Two  Cascaded  Terahertz  Optical  Asymmetric  Demultiplexers  (TOAD), 

In  Proceedings  of  CISS  Conference,  Princeton,  NJ,  USA.  In  press. 

20.  L.  Xu,  B.  C.  Wang,  V.  Baby,  I.  Glesk,  and  P.  R.  Prucnal,  Performance-Improved  All 
Optical  RZ  to  NRZ  Format  Conversion  for  Hybrid  Optical  Networking  Interface,  In 
Proceedings  of  CISS  Conference,  Princeton,  NJ,  USA.  In  press. 

21.  B.  Wang,  L.  Xu,  V.  Baby,  I.  Glesk,  and  P.  Prucnal,  Wavelength-Tunable  Optical 
Spectral  Bistability  in  a  Ring  Laser  Using  a  Semiconductor  Optical  Amplifier,  IEEE 
LEOS'01,  paper  ThLl,  pp.  748-749,  Vol  2,  Hyatt  Regency,  La  Jolla,  San  Diego,  CA,  12 
13  November  2001. 

22.  L.  Xu,  B.  C.  Wang,  V.  Baby,  I.  Glesk,  P.  R.  Prucnal,  M.  Yao,  and  J.  Zhang,  All-Optical 
Clock-Division  Operation  Using  a  Mode-Locked  Fiber  Laser  Based  on  the  Slow  Carrier 
Recovery  Rate  of  SOA,  IEEE  LEOS'01,  paper  MK4,  pp.  87-88,  Vol  1,  Hyatt  Regency, 

La  Jolla,  San  Diego,  CA,  12-13  November  2001. 

23.  L.  Xu,  B.  C.  Wang,  V.  Baby,  I.  Glesk,  and  P.  R  Prucnal,  Widely  Tunable  Fiber  Ring 
Laser  with  EDFA/SOA,  IEEE  LEOS'01,  paper  WC4,  pp.  421-422,  Vol  2,  Hyatt  Regency,  La 
Jolla,  San  Diego,  CA,  12-13  November  2001. 

24.  I.  Glesk,  B.  C.  Wang,  L.  Xu,  D.  Zhou,  R.  J.  Runser,  P.  R.  Prucnal,  Bit  Pattern 
Dependent  Effect  Reduction  in  All-Optical  3R  Regeneration  in  Conference  Proceedings'  3rd 
Electronic  Circuits  and  Systems  Conference,  Bratislava,  Slovakia,  September  5-7, 

2001. 
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25.  B.  C.  Wang,  L.  Xu,  D.  Zhou,  R.  J.  Runser,  I.  Glesk,  and  P.  R.  Prucnal,  Pattern 
Dependent  Effect  Reduction  in  3R  Regeneration  Using  a  Three  Wavelength  Sagnac 
Interferometer  in  Photonics  in  Switching,  OSA  Technical  Digest  (Optical  Society  of 
America,  Washington  DC,  2001,  pp.  49/PThC5-l  -  51/PThC5-3. 

26.  S.  Sears,  Z.  Chen,  H.  Martin,  D.  N.  Christodoulides,  and  M.  Segev,  Clustering  in 
Weakly  Correlated  Wavefronts,  paper  TuQ6,  OSA  Annual  Meeting,  Long  Beach,  California, 
October  14-18,  2001. 

27.  M.  I.  Carvalho,  A.  G.  Grandpierre,  D.  N.  Christodoulides,  and  M.  Segev,  Optical 
Spatial  Shock  Waves  in  Photorefractive  Media,  paper  QTh07,  CLEO/QELS  2001, 

Baltimore,  Maryland,  May  6-11,  2001. 

28.  E.  D.  Eugenieva,  D.  N.  Christodoulides,  and  M.  Segev,  Dynamics  of  Two-Dimensional 
Incoherent  Beams  Resting  on  a  Broad  Intensity  Background,  paper  CThKl,  CLEO/QELS 
2001,  Baltimore,  Maryland,  May  6-11,  2001. 

29.  E.  D.  Eugenieva  and  D.  N.  Christodoulides,  Self-Deflection  of  Partially  Spatially 
Coherent  Beams  and  Solitons  in  Photorefractive  Crystals,  paper  CThK3,  CLEO/QELS  2001, 
Baltimore,  Maryland,  May  6-11,  2001. 

30.  Z.  H.  Musslimani,  M.  Soljacic,  M.  Segev,  and  D.  N.  Christodoulides,  Delayed-Action 
Interaction  Between  Two-Dimensional  Composite  Solitons,  paper  Qth04,  CLEO/QELS 
2001,  Baltimore,  Maryland,  May  6-11,  2001. 

31.  C.  Anastassiou,  M.  Soljacic,  M.  Segev,  E.  D.  Eugenieva,  D.  N.  Christodoulides,  D. 

Kip,  Z.  H.  Musslimani,  and  J.  P.  Torres,  Eliminating  the  Transverse  Instabilities  of  Kerr 
Solitons,  paper  QThK5,  CLEO/QELS  2001,  Baltimore,  Maryland,  May  6-11,  2001. 

32.  D.  Kip,  M.  Soljacic,  M.  Segev,  E.  D.  Eugenieva,  and  D.  N.  Christodoulides, 
Modulation  Instability  and  Pattern  Formation  of  Incoherent  Light,  paper  QWD2, 
CLEO/QELS  2001,  Baltimore,  Maryland,  May  6-11,  2001. 

33.  R.  R.  Malendevich,  Fission  of  Narrow  High  Power  Beams  into  Multiple  Solitons  in  a 
Second-Order  Nonlinear  Slab  Waveguide,  QELS  2001  Technical  Digest,  Baltimore,  QThK6. 

34.  L.  Jankovich,  Observation  of  Satial  Modulation  Instability  in  a  Slab  Waveguide  with 
Third  Order  Nonlinearity,  OSA  Annual  2001,  Long  Beach,  October  2001,  TuD3. 

35.  R.  Schiek,  H.  Fang,  and  G.  I.  Stegeman,  Excitation  of  Quadratic  Solitons,  Topical 
Meeting  on  Nonlinear  Guided  Wave  Phenomena,  Clearwater  Florida,  March  2001, 
postdeadline  paper. 

36.  A.  G.  Grandpierre,  D.  N.  Christodoulides,  M.  I.  Carvalho,  and  M.  Segev,  Optical 
Spatial  Shock  Waves  in  Photorefractive  Media,  paper  MC26,  Nonlinear  Guided  Waves  and 
their  Applications  2001,  Clearwater,  Florida,  March  25-28,  2001. 
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37.  E.  D.  Eugenieva,  D.  N.  Christodoulides,  and  M.  Segev,  Dynamics  of  Two-Dimensional 
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not  always  done  because  of  cost  considerations.  If  a  potential  licensee  is  identified  that  is 
willing  to  assume  the  prosecution  costs  of  filing  a  full  patent  application,  then  a  patent 
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